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ABSTRACT 


|  Navy  Operational  Regional  Atmospheric  Prediction  System 

!  (NORAPS)  analyses  and  forecasts  with  80  km  resolution  are 

used  to  investigate  the  rapid  development  of  an  extratrop- 
ical  cyclone  during  28-29  March  1984  over  the  southeastern 
United  States.  The  lateral  transport  of  vorticity  strongly 
contributes  to  the  spin-up  of  the  low-level  vortex. 
Diabatic  heating  and  moistening  rates  are  estimated  within 
500  km  of  the  storm  from  quasi-Lagrangian  heat  and  moisture 
budgets  and  compared  with  the  same  rates  internally 
predicted  within  the  NORAPS  model.  The  diagnosed  maximum 
diabatic  heating  and  moistening  from  analyses  agree  within 
100  mb  in  elevation  and  6  h  in  time  with  the  model-predicted 
values.  The  model -predicted  diabatic  heating  and  moistening 
rates  are  50%  to  100%  greater  than  diagnosed  budget  calcu¬ 
lations  due  to  an  incorrect  specification  of  the  moistening 
and  heating  rates  in  NORAPS.  Compared  to  more  intense  cases 
of  explosive  maritime  cyclogenesis,  the  higher  level  and 
lower  magnitude  of  maximum  diabatic  heating  is  consistent 
with  the  weaker  deepening  rate  and  continental  nature  of 
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Extratropical  maritime  cyclones  pose  extremely  hazardous 
navigational  and  seaworthiness  problems  for  Naval  operations 
and  commercial  shipping.  To  ensure  safe  maritime  operations 
during  heavy  weather  conditions,  advance  planning  is 
required  to  set  the  heavy  weather  bill,  which  includes 
ballasting  the  ship,  respotting  aircraft  to  the  hangar  deck, 
and  securing  every  compartment  for  heavy  weather.  An  accu¬ 
rate  prediction  of  storm  location  and  intensity  is  needed 
36-48  hours  in  advance  to  allow  the  operational  commander 
time  to  manuever  the  task  force  to  avoid  the  predicted 
region  of  high  winds  and  hazardous  seas. 

While  numerical  weather  prediction  models  have  had 
considerable  success  in  the  prediction  of  many  synoptic- 
scale  systems,  they  perform  rather  poorly  in  predicting  the 
explosive  development  of  maritime  cyclones  (Sanders  and 
Gyakum,  1980).  Explosive  cyclogenesis  is  defined  as  the 
decrease  of  24  mb  per  24  hours  in  the  central  sea-level 
pressure  (SLP)  before  adjusting  to  a  latitude,  <p,  by  multi¬ 
plying  by  sin  <p/sin  60°.  The  highest  frequency  of  explo¬ 
sive  maritime  cyclogenesis  cases  in  the  Northern  Hemisphere 
occurs  over  the  western  ocean  boundary  currents  (Gulf  Stream 
and  Kuroshio)  immediately  east  of  continents  (Sanders  and 
Gyakum,  1980),  where  strong  sea-surface  temperature  (  SST) 
gradients  enhance  the  low-level  baroclinity.  The  Navy 
Operational  Regional  Atmospheric  Prediction  System  (NORAPS) 
has  shown  a  better  success  rate  than  the  Navy  Global 
Atmospheric  Prediction  System  (NOGAPS)  in  predicting  the 
location  and  intensity  of  the  cyclone  center  SLP  for  a  given 
set  of  explosive  maritime  cyclogenesis  cases  (  R.  Hodur, 
NEPRF,  personal  communication).  The  reason  for  this  could 
be  that  the  higher  spatial  and  vertical  resolution  of  the 


regional  model  is  able  to  better  resolve  the  physical 
processes  occurring  within  the  storm  environment.  However, 
a  need  still  exists  for  the  improved  prediction  of  these 
explosive  cyclogenesis  events. 

Maritime  cyclones  are  also  very  important  for  air-sea 
interaction  processes.  The  strong  winds  around  these 
cyclones  may  cause  a  deepening  of  the  oceanic  mixed  layer 
due  to  the  downward  transfer  of  momentum  from  the  atmosphere 
to  the  ocean.  This  fact  has  important  implications  for 
acoustic  propagation  in  the  upper  ocean,  which  anti¬ 
submarine  warfare  (ASW)  and  submarine  commanders  may  use  to 
tactically  exploit  their  situation.  The  strong  cold  advec- 
tion  to  the  rear  of  the  maritime  cyclone  can  produce  posi¬ 
tive  temperature  gradients  ( temperature  increasing  with 
depth)  in  the  upper  ocean  and  lead  to  strong  surface  ducting 
of  acoustic  energy  which  also  can  be  exploited  tactically. 
In  the  future,  regional  oceanographic  models  may  need  to  be 
coupled  to  regional  atmospheric  models  to  predict  correctly 
the  ocean  thermal  structure  response  to  atmospheric  forcing. 

This  thesis  is  a  study  of  an  extratropical  cyclone  which 
formed  over  the  southern  United  States  on  00  GMT  28  March 
1984  and  moved  rapidly  off  the  eastern  coast  of  the  United 
States  by  12  GMT  29  March  1984.  Although  this  storm  is  not 
purely  maritime  in  nature  and  barely  meets  the  criterion  for 
explosive  cyclogenesis,  it  can  serve  as  a  prototype  for 
future  maritime  cases  that  were  observed  in  the  Genesis  of 
Atlantic  Lows  (GALE)  field  experiment  during  1986.  GALE, 
which  was  jointly  sponsored  by  the  Office  of  Naval  Research 
(ONR),  National  Science  Foundation  (NSF)  and  National 
Aeronautics  and  Space  Administration  (NASA),  is  expected  to 
provide  improved  data  sets  over  land  and  the  ocean  for  study 
of  maritime  cyclogenesis.  The  relatively  dense  upper-air 
network  over  land  in  GALE  will  provide  observations  of  the 
horizontal  and  vertical  structure  of  the  atmosphere.  This 


pre-GALE  cyclone  case  was  previously  studied  by  Toll  (1986), 
who  used  a  linear  model  to  investigate  the  instability 
mechanism  responsible  for  the  cyclone  growth.  He  concluded 
that  a  mixed  barotropic-baroclinic  instability  was  respon¬ 
sible  for  the  growth  of  a  mesoscale  perturbation  imbedded 
within  the  larger  scale  cyclonic  circulation. 

The  emphasis  of  this  thesis  will  be  on  the  diabatic 
aspects  involved  in  the  cyclone  growth  and,  in  particular, 
the  role  of  thermodynamic  and  moisture  effects.  Diabatic 
processes  can  provide  an  additional  energy  source  on  the 
smaller  scales  that  may  account  for  the  rapid  intensifica¬ 
tion  of  the  cyclone.  The  NORAPS  analyses  and  forecasts 
constitute  the  data  set  for  various  budget  calculations.  A 
00  GMT  28  March  1984  base  time  is  used  for  the  model  run. 
The  NORAPS  forecast  fields  are  available  at  six  hourly 
intervals  from  00  GMT  28  March  until  12  GMT  29  March  1984. 
The  NORAPS  analyses  are  only  available  every  12  hours  during 
this  period. 

Quasi -Lagrangian  Diagnostics  ( QLD ) ,  which  was  originally 
developed  by  Johnson  and  Downey  (1975)  and  applied  to  numer¬ 
ical  forecasts  of  extratropical  cyclone  cases  by  Wash 
(1978),  are  applied  to  examine  heat,  moisture,  mass  and 
vorticity  budgets.  Several  previous  theses  (Calland,  1983; 
Cook,  1983)  have  used  the  QLD  technique  to  explore  vorticity 
and  mass  budgets  of  various  maritime  cyclones.  Bosse  ( 1984) 
studied  the  role  of  diabatic  effects  on  explosive  maritime 
cyclogenesis  using  NOGAPS,  which  has  a  coarser  spatial  and 
vertical  resolution  than  NORAPS.  These  budget  studies 
investigate  the  mean  properties  of  the  cyclone  both 
spatially  and  temporally,  and  determine  the  relative  contri¬ 
butions  of  the  terms  at  various  stages  of  cyclone  growth. 
Inherent  in  this  technique  is  a  great  amount  of  areal, 
vertical  and  time  averaging  over  the  budget  volume  which  is 
centered  on  the  cyclone.  Because  subdivisions  are  not  made 


within  the  cyclone,  the  results  of  the  budget  study  will  not 
reveal  processes  within  sub-areas  of  the  cyclone. 

The  QLD  technique  uses  an  isobaric  spherical  coordinate 
system  which  is  translated  with  the  storm.  Since  the  radius 
of  the  storm  is  small  compared  to  the  radius  of  the  earth, 
the  budget  volume  can  be  approximated  by  a  cylinder. 
Horizontal  advection  associated  with  storm  translation  is 
isolated  from  the  advection  in  the  moving  coordinate  system, 
so  the  divergence  of  the  transport  (flux)  can  be  associated 
with  cyclone  development  processes.  Vertical  distributions, 
lateral  exchanges  and  sources  and  sinks  of  cyclone  proper¬ 
ties  resulting  from  purely  developmental  processes  are  then 
analyzed. 

This  thesis  is  part  of  a  larger  investigation  into  the 
nature  and  physical  processes  of  maritime  extratropical 
cyclogenesis,  which  has  an  overall  objective  of  improving 
numerical  weather  prediction  over  the  ocean.  The  objectives 
of  this  thesis  are: 

•  Investigate  the  various  properties  of  the  NORAPS  model 
and  identify  areas  where  the  model  does  not  realisti¬ 
cally  represent  the  atmosphere; 

•  Document  the  mean  thermal  and  moisture  structure  of  an 
explosive  cyclogenesis  case  using  QLD  for  both  the 
forecasts  and  analyses; 

•  Assess  the  relative  contribution  of  the  terms  of  the 
thermodynamic  and  moisture  equations  at  various  stages 
during  cyclone  growth,  and  the  horizontal  and  vertical 
distributions  of  those  terms; 

•  Compare  the  moisture  budget  estimates  of  the  diabatic 
heating  rates  with  the  heat  budget  for  the  forecast 
case;  and 

•  Compute  mass  and  vorticity  budgets  during  the  cyclone's 
evolution  to  determine  the  vertical  structure  and 
dynamical  contributions  to  rapid  cyclogenesis. 

A  brief  survey  of  the  literature  on  the  role  of  latent 
heat  release  as  it  affects  extratropical  cyclones  is  given 
in  Chapter  II.  A  synoptic  overview  of  the  NORAPS  analyses 
and  forecasts  in  conjunction  with  other  Numerical  Weather 
Prediction  ( NWP )  models  is  presented  in  Chapter  III.  Mass 
and  vorticity  budget  results  are  discussed  in  Chapter  IV. 


The  central  focus  of  this  thesis  is  contained  in  Chapter  V 
where  the  heat  and  moisture  budget  results  are  discussed. 
An  important  discussion  is  found  in  Appendix  C  which  summa¬ 
rizes  revised  heat  and  moisture  budget  results  after  incor¬ 
porating  corrections  to  a  test  version  of  the  NORAPS  model. 
Conclusions  and  recommendations  for  further  study  are 
outlined  in  Chapter  VI. 


II. 


LITERATURE  SURVEY 

A.  GENERAL 

Thermal  energy  and  differential  vorticity  advection  are 
important  constituents  in  the  growth  of  an  extratropical 
cyclone.  Concentration  of  thermal  advection  forces  greater 
vertical  velocities  and  contributes  to  storm  intensifica¬ 
tion.  This  is  primarily  due  to  a  conservative  adiabatic 
type  process.  However,  the  other  component  of  thermal 
energy  is  the  diabatic  heating  term,  which  if  concentrated, 
can  promote  storm  intensification.  As  diabatic  heating  is 
not  directly  measured  in  the  atmosphere,  heat  and  moisture 
budgets  are  a  technique  of  inferring  its  contribution. 
Diabatic  effects  include  surface  sensible  heat  fluxes, 
latent  heat  release  and  radiation.  How  these  effects  are 
parameterized  in  a  particular  model  is  critical  to  the 
prediction  of  cyclone  development.  These  physical  effects 
must  be  studied  to  gain  a  clearer  understanding  of  the  mari¬ 
time  extratropical  cyclone.  Maritime  cyclones  generally 
have  a  larger  moisture  source  at  the  lower  boundary  which 
contributes  to  earlier  or  more  rapid  intensification  over 
the  ocean  than  over  land.  Two  different  methodologies  may 
be  used  to  study  the  impact  of  moisture  and  latent  heating 
on  cyclone  growth:  sensitivity  studies  using  numerical 

models  and  diagnostic  studies  using  real  or  model-simulated 
data. 

B.  NUMERICAL  STUDIES 

Theories  prior  to  1960  regarding  latent  heat  release 
were  mainly  qualitative  in  nature.  Danard  ( 1964)  was  one  of 
the  first  to  estimate  the  effect  of  latent  heat  release  on 
vertical  velocity  through  horizontal  variations  in  the 
static  stability.  The  results  indicated  the  latent  heat 
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release  served  to  increase  the  vertical  velocity  of  rising 
air  parcels.  He  conjectured  that  latent  heat  release  was  a 
mechanism  for  cyclone  intensification  rather  than  initia¬ 
tion,  since  a  pre-existing  disturbance  had  to  be  present  for 
the  deep  convection  to  occur.  Danard  ( 1966)  later  modeled 
the  vertical  distribution  of  latent  heat  release  as  a  para¬ 
bolic  function  with  a  maximum  in  the  middle  of  the  tropo¬ 
sphere.  Anthes  et  al.  (1983)  and  Gyakum  (1983b)  suggest 
that  a  lower  tropospheric  maximum  in  the  vertical  heating 
profile  is  more  favorable  for  cyclone  intensification. 
Qualitatively,  this  is  consistent  since  the  lower  maximum  in 
vertical  velocity  equates  to  a  stronger  horizontal  conver¬ 
gence  in  the  lower  levels,  which  translates  into  a  greater 
rate  of  low-level  vorticity  generation. 

Gall  (1976)  investigated  the  effects  of  latent  heat 
release  in  growing  baroclinic  waves  with  dry  and  moist 
versions  of  a  general  circulation  model  (GCM).  He  found  the 
growth  rate  of  the  cyclone- scale  wave  (wavenumber  15)  was 
approximately  doubled  in  the  moist  experiment  compared  to 
the  growth  rate  in  the  dry  experiment.  Furthermore,  the 
structure  of  wavenumber  15  was  quite  different  in  the  two 
experiments.  At  the  time  wavenumber  15  reached  its  maximum 
development,  the  moist  model  had  maxima  of  kinetic  energy  at 
500  mb  and  at  the  earth's  surface,  while  the  dry  model  had  a 
single  maximum  at  the  earth's  surface.  The  difference  in 
the  kinetic  energy  spectrum  between  the  two  atmospheres  must 
be  explained  by  the  release  of  latent  heat.  Gall  (1976) 
attributed  this  difference  to  the  moist  convective  adjust¬ 
ment  by  which  the  temperatures  at  the  middle  and  upper- 
levels  are  increased  relative  to  those  at  the  earth' s 
surface.  He  also  found  a  difference  in  the  500  mb  tempera¬ 
ture  perturbations  east  and  west  of  the  most  intense  surface 
lows.  In  the  moist  experiment,  the  magnitude  of  the  500  mb 
warm  temperature  perturbation  east  of  these  lows  exceeded 


the  magnitude  of  the  cold  temperature  perturbation  west  of 
these  lows  by  2  to  5°K,  whereas  the  perturbations  on  either 
side  of  these  lows  were  approximately  the  same  in  the  dry 
experiment. 

Mak  (1982)  discussed  moist  quasi-geostrophic  (Q-G)  baro- 
clinic  instability  by  generalizing  the  analytic  theory  of 
Q-G  dynamics  in  terms  of  condensational  heating.  Since  the 
time  and  length  scales  of  the  precipitation  pattern  associ¬ 
ated  with  a  typical  extratropical  cyclone  were  comparable  to 
those  of  the  primary  circulation,  there  was  likely  a  close 
direct  feedback  between  the  heating  and  circulation.  This 
formulation  was  somewhat  similar  to  the  Conditional 
Instability  of  the  Second  Kind  (CISK)  theory  proposed 
jointly  by  Charney  and  Eliassen  (1964)  and  Ooyama  (1964). 
Mak' s  premise  was  the  moisture  supply  for  the  condensational 
process  in  a  baroclinic  wave  was  sustained  by  the  low-level 
convergence  field  of  the  wave.  The  condensational  heating 
was  parameterized  in  terms  of  the  vorticity  field  of  the 
disturbance  rather  than  the  vertical  velocity.  This  parame¬ 


terization  was  more  effective  for  low  frequency  disturbances 
or  those  with  the  longer  wavelengths.  Mak's  results 
revealed  that  as  the  heating  intensity  parameter  was 
increased,  the  growth  rate  of  the  unstable  wave  signifi¬ 
cantly  increased,  the  wavelength  significantly  decreased  and 
the  phase  speed  increased.  His  results  supported  the 
general  notion  that  the  baroclinic  forcing  in  a  disturbance 
can  organize  the  condensational  heating  on  a  scale  compa¬ 
rable  to  the  wave. 

Chang  et  al.  (1982)  have  simulated  the  cyclone  system 
with  and  without  latent  heating  using  a  regional  fine  mesh 
(  140  km)  model  encompassing  the  United  States.  In  the  dry 
simulation,  the  model  failed  to  predict  the  formation  of  a 
closed  circulation  throughout  the  depth  of  the  troposphere 
and  a  pronounced  northwest- southeast  horizontal  tilt  of  the 
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upper-level  trough.  The  moist  simulation  produced  both  of 
these  features  which  were  in  good  agreement  with  NMC  anal¬ 
yses.  The  significance  of  the  latent  heat  release  was 
illustrated  by  subtracting  the  predicted  dry  values  from  the 
corresponding  moist  values.  At  500  and  300  mb,  the  largest 
positive  temperature  differences  were  confined  to  the  areas 
of  large  precipitation,  which  indicates  that  latent  heating 
was  the  main  factor  contributing  to  this  positive  tempera¬ 
ture  difference.  However,  the  moist  model  had  lower  700  mb 
temperatures  than  the  dry  model.  Their  results  indicated 
that  latent  heat  release  stabilized  the  middle  and  upper 
troposphere  and  reduced  the  large-scale  horizontal  tempera¬ 
ture  gradient  or  baroclinity.  They  also  postulate  that  the 
increased  circulation  in  the  lower  troposphere  caused  the 
formation  of  a  low-level  jet. 

Chang  et  al.  (1984)  also  investigated  the  latent  heat 
induced  energy  transformations  during  cyclogenesis. 
Although  condensational  processes  were  not  efficient  in 
increasing  the  total  potential  energy  of  the  model  atmos¬ 
phere,  latent  heat  acted  as  a  catalyst  to  enhance  the 
conversion  of  potential  to  kinetic  energy  within  the  mid¬ 
latitude  cyclone.  The  most  significant  response  to  heating 
appeared  in  the  lower  troposphere.  Latent  heating  changed 
the  configuration  of  the  motion  field  from  a  large-scale 
pattern  with  pronounced  vertical  shear  into  a  small-scale 
feature  with  marked  horizontal  shear  which  served  to  reduce 
the  sub-grid  scale  dissipation  below  500  mb.  Additionally, 
the  enhancement  of  ageostrophic  generation  of  kinetic  energy 
provided  an  important  source  of  energy  for  the  maintenance 
of  the  lower  tropospheric  circulation. 

Sardie  and  Warner  (1983)  included  the  effects  of  latent 
heating  and  baroclinity  in  an  analytic  three-layer  Q-G  model 
to  investigate  the  mechanism  of  polar  low  development.  This 
mesoscale  phenomenon  has  wavelengths  of  500-1000  km,  has  a 
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range  of  vertical  scales  and  develops  within  the  cold  air 
behind  a  cold  front.  They  integrated  their  model  for  seven 
real-data  cases  using  three  modes,  including:  a  dry  baro- 
clinic  case  (no  latent  heading);  a  moist  baroclinic  case  in 
which  diabatic  heating  is  partitioned  between  the  convective 
and  non-convecti ve  parts;  and  the  case  in  which  latent 
heating  is  solely  caused  by  convective  clouds,  referred  to 
as  CI3K.  Their  results  indicate  Atlantic  polar  lows  develop 
due  to  a  T 1 .7  >  type  mechanism  shallow  baroclinity,  while 


moist  bar 


is  the  iominaut  mechanism  for  Pacific 


pclar  l.vs.  T:.e:r  results  on  the  Atlantic  polar  lows  are 
■:c:is  i  ste;;'  w :  *. :.  Pasmussen  (  17"?),  who  argues  that  some  type 
of  CIS.-'  met  .nanism  must  be  operating  for  these  lows  which 
have  the  smallest  horizontal  scales. 

Anth.es  et  a..  (  1983)  performed  a  series  of  numerical 
simulations  using  the  Queen  Elizabeth  II  storm.  The  results 
from  the  latent  heating  investigation  indicated  the  evolu¬ 
tion  of  the  model  storm  was  not  significantly  affected  by 
latent  heating  in  the  early  stages.  As  the  storm  intensi¬ 
fied  and  vertical  motions  became  stronger,  latent  heating 
played  a  greater  role  in  enhancing  the  development.  Gyakum 
(1983)  also  demonstrated  that  quasi-geostrophic  dynamics 
could  not  account  for  the  observed  intensity  of  the  cyclone 
during  the  explosive  stage.  He  reasoned  that  the  diabatic 
heating  in  this  case  was  responsible  for  the  extreme  thick¬ 
ness  changes  in  the  lower  part  of  the  column.  He  calculated 
the  residual  warming  of  the  column  over  the  surface  low 
center  to  be  16°C  in  a  12-hour  period. 

C.  BUDGET  AND  DIAGNOSTIC  STUDIES 

The  purpose  of  budget  studies  is  to  diagnose  terms  such 
as  diabatic  heating  and  precipitation  by  properly  accounting 
for  all  the  sources  and  sinks  through  the  calculation  of 
terms  in  the  thermodynamic  and  moisture  conservation  equa¬ 
tions.  Kuo  and  Anthes  (1984)  state  that  budget  calculations 


are  subject  to  errors  due  to  sampling,  measurement,  finite 
differencing  and  objective  analysis.  Errors  associated 
with  wind  measurements  can  introduce  substantial  errors  in 
the  calculation  of  divergence  and  vertical  velocity,  which 
are  two  of  the  most  crucial  terms  in  the  budget  calcula¬ 
tions.  Another  factor  which  can  introduce  errors  is  spatial 
and  temporal  interpolation  performed  during  the  objective 
analysis.  Therefore,  a  detailed  error  analysis  is  of 
particular  importance  to  budget  studies,  since  the  residual 
of  the  calculation  is  influenced  by  all  of  the  errors 
mentioned  above.  The  residual  also  includes  real,  sub-grid 
scale  physical  effects,  which  can  be  interpreted  only  after 
the  contribution  of  errors  is  isolated. 

Kuo  and  Anthes  (  1984)  note  that  the  mean  divergence 
removal  technique  (O'  Brien,  1970)  generally  has  a  stronger 
influence  on  the  heat  budget  than  the  moisture  budget. 
Because  the  moisture  decreases  rapidly  with  height,  erro¬ 
neous  vertical  motions  at  upper  levels  cause  smaller  errors 
in  the  moisture  budget.  In  their  study,  they  found  area- 
averaged  errors  in  the  heat  and  moisture  budgets  of  about 
5°C  and  2  g/kg,  respectively,  for  a  temporal  scale  of  6 
hours  and  a  spatial  scale  of  550  x  550  km.  The  vertically- 
integrated  changes  due  to  radiation,  sensible  heat  flux  and 
moisture  changes  due  to  evaporation  are  small  compared  to 
the  inherent  root  mean  square  (rms)  errors  in  the  budget 
calculation.  In  general,  the  above  errors  are  small 
compared  to  those  due  to  latent  heating  during  the  mature 
stage  of  convection,  which  are  typically  20-30°C/day  and  10 
g/ kg/da y. 

Smith  et  al.  ( 1984)  established  the  significance  of 
latent  heat  release  in  developing  extratropical  cyclones  by 
separating  the  vertical  motion  forced  by  latent  heat  release 
from  the  total  vertical  motion.  The  total  vertical  motion 
was  calculated  using  the  kinematic  method,  while  the 


vertical  motion  due  to  latent  heat  release  was  determined 
from  the  omega  equation.  The  ratio  of  the  vertical  motion 
due  to  latent  heat  to  the  total  vertical  motion  was  then 
calculated.  In  their  case  study,  a  mid-period  maximum  of 
90%  below  900  mb  suggested  the  increased  latent  heat  release 
had  a  particularly  strong  influence  on  below-cloud  vertical 
motions.  Additionally,  a  comparison  was  made  of  the 
vertically-integrated  convective  and  stable  latent  heat 
release.  The  total  latent  heat  release  was  dominated  by  the 
convective  components,  even  though  this  was  a  mid-winter 
storm  over  the  north-central  United  States.  A  key  finding 
was  that  the  maximum  rate  of  development  occurred  after  the 
time  of  maximum  latent  heating. 

Liou  and  Elsberry  (  1985)  performed  a  QLD  study  of  an 
explosive  maritime  cyclogenesis  case  over  the  northwest 
Pacific  Ocean  using  a  research  version  of  the  Univerisity  of 
California  at  Los  Angeles  (UCLA)  GCM.  Their  results 
revealed  area-averaged  heating  rates  of  25-30°C/day  with  the 
maximum  level  at  600-700  mb  and  that  diabatic  heating  was 
caused  by  latent  heat  release  from  stable  condensation  and 
middle-level  convection.  Additionally,  they  found  that  the 
sea-level  pressure  (SLP)  deepening  rate  was  highly  corre¬ 
lated  with  the  diabatic  heating  rate. 

A  diagnostic  moisture  analysis  of  a  meso-|5  scale  thun¬ 
derstorm  environment  was  performed  by  Fuelberg  et  al. 
(1986).  A  highly  refined  rawindsonde  network  of  75  km 
spacing  (over  Oklahoma)  with  3  hourly  soundings  was  used.  A 
stationary  15  x  13  budget  grid  was  formed  with  25  km  hori¬ 
zontal  spacing.  Data  were  available  every  50  mb  from  900  to 
150  mb  with  the  exception  of  specific  humidity  which  was  not 
reported  above  350  mb.  Coarser  resolution  (synoptic)  data 
were  available  from  the  National  Weather  Service  (NWS) 
rawinsonde  network  which  allowed  a  separate  budget  calcula¬ 
tion  for  comparison  with  the  finer  resolution  data.  Their 
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results  showed  that  horizontal  moisture  convergence  was  the 
most  important  term  in  the  budget,  especially  in  the  surface 
to  850  mb  layer.  Most  of  the  horizontal  flux  convergence  of 
moisture  was  attributed  to  the  velocity  convergence,  with 
moisture  advection  being  of  secondary  importance.  The  hori¬ 
zontal  moisture  convergence  increased  an  order  of  magnitude 
when  the  storms  were  occurring.  The  residuals  indicated 
vapor  accumulation  prior  to  the  storm  activity,  but  reversed 
to  condensation/precipitation  once  the  storms  began.  The 
residuals  also  increased  an  order  of  magnitude  as  in  the 
horizontal  flux  term.  A  comparison  with  the  budget  based  on 
the  synoptic  data  showed  the  magnitude  of  terms  in  the 
mesoscale  budget  approached  an  order  of  magnitude  greater 
than  the  synoptic-scale  values  near  the  time  of  the  storms. 


A.  GENERAL 

A  discussion  of  various  model  forecasts  and  verifying 
analyses  is  presented  to  document  the  synoptic  evolution  of 
the  28-29  March  1984  storm.  This  storm  is  frequently 
referred  to  as  the  Carolinas  storm  as  a  particularly  severe 
outbreak  of  tornadoes  caused  widespread  damage  in  North  and 
South  Carolina.  Various  meteorological  aspects  will  be 
highlighted  and  physical  mechanisms  presented  to  explain  the 
growth  of  the  cyclone.  Significant  differences  between  the 
model  forecast  and  verifying  analyses  will  be  mentioned. 
For  convenience  of  comparison,  the  figures  are  grouped  at 
the  end  of  the  chapter.  In  this  study,  the  NORAPS  surface 
analyses  are  used  as  the  basis  for  verifying  the  corre¬ 
sponding  NORAPS  forecast.  The  NORAPS  surface  analysis 
procedure  includes  using  an  89  x  89  hemispheric  grid  (268 
km)  to  interpolate  to  the  higher  resolution  (80  km)  NORAPS 
grid.  A  regional  update  cycle  uses  the  12  h  NORAPS  forecast 
fields  to  provide  the  "first  guess"  for  the  NORAPS  forecast. 
The  National  Meteorological  Center  (NMC)  analyses,  in 
particular,  provide  a  different  interpretation  since  they 
are  hand-drawn  and  subjective  in  nature  as  opposed  to  NORAPS 
objective  analyses.  These  NMC  analyses  also  depict  much 
more  of  the  mesoscale  detail  within  the  larger  scale 
cyclonic  circulation.  Mesoscale  low  centers,  which  have  a 
horizontal  scale  of  100-200  km,  are  very  evident  in  the 
series  of  NMC  analyses  for  this  storm.  Toll  (1986)  used 
rawindsonde  observations  and  a  NORAPS  analysis  at  00  GMT  29 
March  1984  to  construct  cross-sections  for  a  linear 
stability  analysis  to  investigate  the  instability  mechanisms 
associated  with  cyclone  growth. 
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The  tracks  of  the  storm  for  both  the  NORAPS  forecast 
(every  6  h)  and  analyses  (every  12  h)  from  00  GMT  28  March 
1984  until  12  GMT  29  March  1984  are  presented  in  Fig.  3. la. 
The  model-predicted  track  of  the  cyclone  agrees  well  with 
the  analyzed  position  until  00  GMT  29  March  1984.  During 
the  following  12  h,  the  model-predicted  track  departs 
significantly  from  the  analyzed  track  of  the  cyclone  to  the 
south.  The  NORAPS  predicted  cyclone  intensity  (Fig.  3.1b) 
agrees  quite  well  with  the  analyzed  intensity,  especially  at 
the  later  time  periods. 

B.  NORAPS  ANALYSIS  -  00  GMT  28  MARCH  1984 

The  initial  surface  analysis  (Fig.  3.2a)  at  00  GMT  28 
March  1984  depicts  a  broad  990  mb  center  of  low  pressure 
over  eastern  Texas.  The  corresponding  1000-500  mb  thickness 
pattern  (Fig.  3.2b)  reveals  moderate  warm  advection  over  the 
southeastern  U.  S.  and  strong  cold  advection  across  the 
southwestern  U. S.  The  low-level  warm  advection  patterns  at 
925  and  850  mb  (not  shown)  have  two  distinct  thermal 
gradient  zones:  one  along  the  Gulf  Coast  and  the  other  is 
displaced  much  farther  north  through  the  Ohio  Valley.  This 
lack  of  concentrated  low-level  baroclinity  in  advance  of  the 
storm  is  consistent  with  a  weak  upward  vertical  motion 
pattern.  A  500  mb  pressure  trough  (Fig.  3.2b)  is  located 
along  100°W  and  a  thermal  trough  lags  the  height  trough  by  a 
quarter  of  a  wavelength.  This  pattern  satisfies  a  necessary 
condition  for  baroclinic  instability.  The  broad  elongated 
trough  at  500  mb  indicates  a  fairly  weak  relative  vorticity 
pattern  due  to  curvature  effects  alone.  The  925,  850  and 
700  mb  analyses  (not  shown)  depict  a  short-wave  trough  near 
35°N,  110°W  (west  of  the  500  mb  trough  axis)  that  rotates 
through  the  elongated  trough  at  500mb.  The  wind  analysis  at 
300  mb  (Fig.  3.3)  has  a  60  m/s  jet  streak  in  the  vicinity  of 
western  Louisiana.  The  strong  cyclonic  shear  north  of  this 
jet  streak  contributes  heavily  to  the  positive  relative 


vorticity  at  300  mb  since  the  curvature  component  is  weak  in 
this  region  as  noted  above.  The  surface  low  center  is 
located  in  the  right-rear  quadrant  of  this  jet,  which  is 
favorable  for  surface  deepening  due  to  upper-level 
divergence. 

C.  NWP  PERFORMANCE 

NORAPS  forecasts  and  verifying  analyses  will  be  empha¬ 
sized  in  this  section.  Other  models  and  analyses  will  be 
introduced  during  the  discussion  to  highlight  the  high 
degree  of  variability  found  in  these  products.  Physical 
interpretations  for  the  differences  among  the  various  models 
will  be  given.  It  is  important  to  note  that  a  25-point 
filter  is  applied  eight  times  to  the  NORAPS  forecast  output 
fields  for  display  purposes,  which  may  tend  to  smooth  out 
some  of  the  forecast  mesoscale  structure.  A  technique 
called  Fields  by  Information  Blending  (FIB)  is  performed  to 
produce  the  NORAPS  surface  analyses  using  an  89  x  89  hemis¬ 
pheric  grid.  This  technique  also  results  in  a  smoothed 
analysis. 

The  NORAPS  model  forecast  (Fig.  3.4a)  at  06  GMT  28  March 
1984  fills  the  surface  low  center  4  mb  to  994  mb  as  it  moves 
northeast  to  33°N,  90°W.  Although  a  corresponding  NORAPS 
analysis  is  not  available  to  verify  the  filling  of  the  low 
during  this  6  h  period,  a  NMC  analysis  valid  at  this  time 
does  not  indicate  this  filling  tendency  (taking  into  account 
the  diurnal  variation  in  pressure).  The  low-level  warm 
advection  pattern  (not  shown)  is  rather  diffuse  while  strong 
cold  air  advection  continues  over  Texas.  The  500  mb  fore¬ 
cast  (Fig.  3.4b)  depicts  a  more  elongated  trough  than  in  the 
analysis.  This  is  a  consequence  of  the  short  wave  trough 
discussed  above  having  translated  into  the  western  part  of 
the  trough.  The  predicted  300  mb  wind  maximum  (not  shown ) 
has  weakened  by  about  10  m/s  from  the  analysis.  The  filling 
tendency  predicted  by  the  NORAPS  forecast  may  be  attributed 


to  adjustments  in  the  initial  mass  and  momentum  fields  which 
occur  in  the  first  few  hours  after  the  time  integration  of 
the  model. 

At  12  GMT  28  March  1984,  the  NORAPS  surface  forecast 
(Fig.  3.5a)  depicts  a  990  mb  low  center  near  36°N,  87°W, 
whereas  the  corresponding  analysis  (Fig.  3.5b)  has  a  987  mb 
center  near  36°N,  89°W.  This  3  mb  pressure  difference 
between  the  analysis  and  forecast  is  the  largest  difference 
during  the  forecast.  Both  the  model  and  analyzed  low 
tropospheric  fields  depict  two  distinct  thermal  gradient 
regions  in  advance  of  the  low  center  which  form  into  a 
single  zone  to  the  rear  of  the  low.  The  NORAPS  500  mb  fore¬ 
cast  (Fig.  3.6a)  underestimates  the  intensity  of  the  trough 
by  60  m.  The  300  mb  wind  field  (not  shown)  verifies  well 
with  a  60  m/s  jet  over  southern  Louisiana  in  both  the  fore¬ 
cast  and  analysis. 

In  the  18  GMT  28  March  1984  N0RAPS  forecast  (Fig.  3.7a), 
the  surface  low  deepens  to  987  mb  and  translates  eastward  to 
36°N,  85°W.  Greater  organization  is  seen  in  the  low-level 
thermal  field  (not  shown)  in  advance  of  the  low  which  indi¬ 
cates  the  formation  of  a  warm  front  in  the  lower  tropo¬ 
sphere.  The  500  mb  trough  axis  translates  to  95°W  with  60  m 
height  falls  occurring  at  the  base  of  the  trough  (  F:,.g. 
3.7b).  The  300  mb  wind  field  (Fig.  3.8)  indicates  an  inter¬ 
esting  pattern  of  an  intensifying  jet  streak  ( in  the 
northwesterly  flow)  associated  with  an  upper-level  trough 
downstream  of  the  major  upper-level  trough  being  investi¬ 
gated.  Another  60  m/s  jet  streak  is  associated  with  the 
major  upper-level  trough  over  the  Gulf  of  Mexico.  The 
developing  surface  low  is  in  the  right-rear  and  left-front 
quadrants  of  the  two  jet  streaks  previously  mentioned.  The 
expected  divergence  in  these  quadrants  of  the  jets  can 
enhance  the  surface  development  through  removal  of  mass 


In  the  00  GMT  29  March  1984  NORAPS  forecast  (Fig.  3.9a), 
the  surface  low  deepens  to  983  mb  and  tracks  eastward  to 
81°W.  The  verifying  analysis  (Fig.  3.9b)  shows  excellent 
agreement  with  the  model  with  respect  to  position  and  inten¬ 
sity  (Fig.  3.1).  The  1000-500  mb  thickness  pattern  (Fig. 
3.9a  and  Fig.  3.9b)  indicates  strong  cold  advection  along 
the  Gulf  Coast  and  strong  warm  advection  along  the  eastern 
seaboard  north  of  Cape  Hatteras  in  both  the  forecast  and 
analysis.  The  predicted  500  mb  trough  does  not  verify  well 
as  the  central  value  is  60  m  too  high  and  the  formation  of  a 
closed  low  is  not  predicted  (Fig.  3.10a).  Whereas  a  60  m/s 
jet  streak  is  predicted  (Fig.  3.11a)  along  the  Gulf  Coast, 
the  analyzed  maximum  is  actually  80  m/s  (Fig.  3.11b).  The 
extensive  outbreak  of  severe  weather,  including  22  torna¬ 
does,  in  the  Carolinas  occurs  within  a  6  h  period  centered 
on  00  GMT  29  March  1984,  as  outlined  by  Ferguson  et.  al. 

( 1986). 

Another  model,  the  Mesoscale  Atmospheric  Simulation 
System  (MASS),  developed  by  NASA,  was  used  by  Kocin  et  al. 
(1984)  to  simulate  the  severe  weather  outbreak  across  the 
Carolinas.  This  model's  smaller  grid  size  (50  km)  and 
extensive  physics  package  can  sometimes  capture  the  dynam¬ 
ical  interactions  and  diabatic  processes  which  force  mesos¬ 
cale  features.  The  model  was  initialized  at  12  GMT  28  March 
1984  and  integrated  15  hours  (to  03  GMT  29  March  1984)  to 
study  the  evolution  of  a  mesoscale  low  that  was  involved  in 
triggering  the  severe  weather.  Forecasts  are  available 
every  3  hours.  A  time  series  of  the  MASS  forecast  and  NMC 
analyses  is  presented  in  Fig.  3. 12.  A  significant  feature 
of  this  model  prediction  was  the  formation  of  a  mesoscale 
low  center  between  18-21  GMT  28  March  1984  over  eastern 
Alabama,  which  moved  rapidly  northeast  and  was  associated 
with  the  severe  weather  event  in  the  Carolinas.  Although 
this  model  prediction  was  3  h  too  slow  in  the  prediction  of 
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this  mesoscale  feature  and  lagged  the  analyzed  low  center  by 
250  km,  it  more  accurately  modeled  the  observed  state  than 
NORAPS.  The  possible  reasons  why  MASS  produced  a  more 
accurate  forecast  than  NORAPS  are  that  its  extensive  physics 
package  and  smaller  grid  size  enabled  it  to  capture  some  of 
the  low-level  forcing  mechanisms  such  as  the  surface  heat 
fluxes  and  latent  heat  release  on  the  smaller  scales,  which 
could  not  be  adequately  resolved  in  the  NORAPS  prediction. 
Kocin  et  al.  ( 1984)  attribute  model  deficiencies  in  MASS  to 
a  poor  initial  data  base  at  jet  stream  level  where  several 
critical  wind  obsevations  were  missing. 

The  06  GMT  29  March  1984  NORAPS  forecast  (Fig.  3.13a) 
continues  the  linear  decrease  in  the  central  sea-level  pres¬ 
sure  to  978  mb.  However,  a  marked  departure  in  the  track 
occurs  (Fig.  3.1)  as  the  storm  is  predicted  to  decelerate 
and  move  southeast.  The  southeast  movement  may  be  a  delayed 
response  in  NORAPS  to  some  mesoscale  forcing  which  formed  a 
mesoscale  low  center  south  of  the  major  low  at  21  GMT  28 
March  1984  (Fig.  3.12).  The  model  cannot  adequately  resolve 
the  mesoscale  forcing  mechanism  and  therefore,  decelerates 
the  major  low  center  and  steers  it  southeast.  A  separate 
integration  of  NORAPS  in  which  the  convective  component  of 
the  precipitation  is  eliminated  has  a  more  rapid  movement  of 
the  low  center  to  the  northeast  at  the  corresponding  time. 
Strong  latent  heat  release  in  a  localized  area  may  tend  to 
decelerate  the  storm.  The  500  mb  NORAPS  forecast  (  Fig. 
3.  13b)  depicts  a  closed  5400  m  low  just  to  the  west  of  the 
surface  low.  The  300  mb  predicted  wind  field  (not  shown) 
continues  to  position  the  jet  streak  along  the  Gulf  Coast, 
although  it  has  probably  translated  east  of  the  trough  axis 
based  on  the  severe  weather  outbreak  in  the  Carolinas.  The 
incorrect  prediction  of  the  translation  of  this  jet  streak 
may  have  also  contributed  to  the  erroneous  surface  low 
deceleration  and  movement  to  the  east- southeast. 
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In  the  12  GMT  29  March  1984  NORAPS  forecast  (Fig. 
3.14a),  the  low  center  of  974  mb  verifies  well  with  respect 
to  the  NORAPS  analyzed  intensity  of  973  mb.  However,  the 
forecast  position  is  360  n  mi  southwest  of  the  analyzed 
position.  The  NMC  surface  analysis  at  this  time  depicts  a 
966  mb  surface  low  center  off  the  eastern  coast  of  Maryland, 
which  correlates  well  with  the  N0RAP3  analyzed  position,  but 
is  7  mb  lower  in  pressure.  The  circulation  of  the  storm  is 
quite  extensive  and  covers  the  entire  eastern  half  of  the  U. 
S.  and  western  Atlantic.  The  low-level  thickness  pattern 
(Fig.  3.14b)  indicates  warm  advection  extends  around  the  low 
center  and  into  the  northwest  quadrant.  The  system  is 
becoming  nearly  vertically  stacked  in  both  the  forecast  and 
analysis,  which  indicates  it  has  reached  maximum  intensity 
and  has  converted  most  of  the  available  potential  energy  to 
kinetic  energy.  Although  the  300  mb  analysis  (not  shown) 
has  a  70  m/s  jet  streak  east  of  Georgia,  a  60  m/s  jet  streak 
remains  quasi-stationary  over  the  northern  Gulf  of  Mexico  in 
the  prediction. 

A  single  36  h  Limited  Fine  Mesh  (LFM)  forecast  valid  at 
12  GMT  29  March  1984  is  available  for  comparison  with  the 
NORAPS  model.  A  974  mb  low  center  is  predicted  by  the  LFM 
model  to  be  in  central  Virginia,  which  is  somewhat  better 
than  the  NORAPS  36  h  predicted  location  in  southern  North 
Carolina  (Fig.  3.14a).  The  LFM  forecast  intensity  of  the 
low  center  (  974  mb)  was  the  same  as  the  NORAPS  forecast. 
This  is  8  mb  too  high  according  to  the  NMC  analysis  but  only 
1  mb  too  high  according  to  the  NORAPS  analysis.  The  LFM  500 
mb  forecast  trough  verified  extremely  well  with  a  predicted 
5300  m  closed  low  center,  whereas  a  5280  m  center  was 
analyzed.  The  accumulated  6  h  precipitation  chart  has  large 
"bullseye"  values  in  excess  of  two  inches.  Both  NORAPS  ( see 
Fig.  5.  14)  and  LFM  models  overpredict  precipitation  for 
this  storm. 
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NOGAPS  is  a  global  model  of  relatively  coarse  resolution 
(2.4°  lat  by  3.0°  long)  that  is  designed  primarily  for 
longer-range  weather  prediction  over  the  open  o:ean.  This 
model  is  not  able  to  resolve  the  mesoscale  detail  which  the 
regional  models  can  achieve.  The  NOGAPS  36  h  forecast  valid 
at  12  GMT  29  March  1984  had  a  976  mb  surface  low  center  over 
western  Virginia,  whereas  the  analyzed  position  was  some  500 
n  mi  farther  east.  The  intensity  of  the  predicted  low  was 
10  mb  too  high  compared  to  the  NMC  analysis  and  3  mb  too 
high  compared  to  the  NORAPS  analysis.  It  is  difficult  to 
determine  whether  numerical  truncation  or  physical  errors 
accounted  for  the  slower  predicted  movement  of  the  cyclone 
in  NOGAPS. 

D.  SUMMARY 

The  intensity  of  the  storm  is  predicted  quite  well  by 
the  NORAPS  model  as  verified  by  the  NORAPS  analyses.  The 
track  of  the  low  is  also  predicted  fairly  well  for  the  first 
24  hours  although  a  large  position  error  is  evident  at  36 
hours.  A  sharp  difference  exists  in  the  low  center  intensi¬ 
ties  at  12  GMT  29  March  1984  between  the  NORAPS  and  NMC 
analyses,  which  most  likely  can  be  attributed  to  the 
smoothing  in  the  NORAPS  objective  analysis  scheme  and  the 
coarse  resolution  of  268  km  in  the  surface  analysis.  It  is 
difficult  to  isolate  whether  the  poor  track  prediction 
during  the  last  12  h  was  a  result  of  NORAPS  incorrectly 
responding  to  some  mesoscale  forcing  mechanism  occurring  to 
the  southwest  of  the  major  low  (as  in  the  MASS  model)  or 
whether  it  overresponded  to  convection  occurring  in  the 
vicinity  of  the  low.  The  organization  of  the  multiple  low 
centers  into  a  consolidated  center  by  12  GMT  29  March  1984 
suggests  orographic  forcing  was  responsible  for  a  large 
portion  of  the  mesoscale  detail  while  the  the  storm  was  in 
the  vicinity  of  the  Appalachians. 
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As  m  Fig.  3.2,  except  for  06  GMT  29  March  196 
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IV.  MASS  AND  VORTICITY  BUDGET  RESULTS 
A.  MASS  BUDGET 

The  mass  budget  is  performed  using  the  finite  difference 
form  of  the  continuity  equation  in  pressure  coordinates. 
The  values  of  vertical  velocity,  commonly  referred  to  as 
omega,  are  computed  kinematically  and  interpolated  to  the 
budget  volume.  Horizontal  winds  are  likewise  interpolated 
to  the  budget  volume  with  normal  and  tangential  wind  compo¬ 
nents  being  computed  for  budget  calculations.  A  method 
developed  by  O'Brien  (1970)  is  used  to  adjust  the  horizontal 
and  vertical  winds  to  achieve  mass  balance.  In  this  method, 
the  total  horizontal  flux  for  a  particular  radius  ring  of 
the  budget  volume  is  computed  using  the  line  integral 
method.  These  horizontal  fluxes  are  summed  in  the  vertical 
to  arrive  at  a  vertically-integrated  horizontal  mass  flux 
for  each  radius  ring.  The  vertically-integrated  horizontal 
flux  is  adjusted  to  balance  the  net  vertical  mass  flux 
between  the  top  (100  mb)  and  an  arbitrary  chosen  bottom 
level  of  1000  mb.  In  the  forecast  case,  the  model-predicted 
values  of  omega  are  used  at  the  upper  and  lower  boundaries, 
while  omega  is  set  equal  to  zero  at  the  boundaries  in  the 
analysis  case. 

An  adjusted  divergence  and  normal  wind  are  then  calcu¬ 
lated  at  36  points  around  the  budget  center  at  each  level. 
Based  on  these  two  adjustments,  an  adjusted  vertical 
velocity  is  computed  at  interior  levels.  For  the  NORAPS 
forecast  fields,  a  constant  correction  factor  is  used  to 
adjust  the  horizontal  divergence  at  each  level.  For  the 
analyses  fields,  a  weighting  function  that  linearly 
decreases  with  pressure  is  used  to  correct  the  divergence 
due  to  the  larger  error  in  the  wind  estimates  with  height. 
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Besides  the  above  adjustments  required  for  mass  balance, 
there  can  be  other  sources  of  errors  in  the  budget  analyses. 
Serious  truncation  error  can  be  introduced  when  using  finite 
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differences  and  interpolating  from  sigma  levels  in  the  model 
to  pressure  levels  in  the  budgets.  Also,  errors  can  arise 
when  interpolating  from  the  Lambert  conformal  grid  to  the 
cylindrical  budget  volume. 

B.  MASS  BUDGET  RESULTS 

The  mass  budget  consists  of  a  horizontal  and  vertical 
mass  transport  term.  The  horizontal  mass  transport  term 
represents  convergence  or  divergence  within  a  particular 
layer  (averaged  over  two  levels)  with  no  time  averaging 
involved.  The  vertical  mass  transport  is  inferred  from  the 
vertical  velocity  field,  which  has  been  adjusted  according 
to  O'Brien's  (1970)  method  discussed  earlier.  The  vertical 


velocity  fields  are  averaged  over  two  time  periods,  however, 
no  layer  averaging  is  involved. 

A  comparison  will  be  made  of  the  NORAPS  forecast  hori¬ 
zontal  and  vertical  transport  fields  with  the  analysis 


fields.  The  forecast  fields  are  available  every  6  h  from  06 
GMT  28  March  until  12  GMT  29  March  1984.  The  analysis 
fields  are  available  only  at  12  h  intervals  from  00  GMT  28 
March  to  12  GMT  29  March  1984.  Mass  fluxes  and  vertical 
velocities  in  a  vertical/time  section  for  radius  4°  lat. 
will  be  emphasized  because  this  radius  most  nearly  repre¬ 
sents  the  inner  core  of  the  cyclone  where  heating  and  mois¬ 
tening  tend  to  be  concentrated. 


The  analyzed  horizontal  mass  transport  (Fig.  4.1a)  has  a 
pattern  of  inward  mass  transport  below  700  mb  and  mass 
outflow  or  divergence  above  700  mb.  The  forecast  horizontal 
transport  (Fig.  4.1b)  has  a  similar  two-layer  vertical 
structure.  However,  the  level  of  non-divergence  (LND), 
where  inward  and  outward  mass  transport  changes  sign,  occurs 
at  500  mb.  The  much  lower  LND  in  the  analyzed  case  is 
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typical  of  most  explosive  maritime  cyclones  which  have  been 
found  to  have  a  shallow  vertical  structure  ( Sanders  and 
Gyakum,  1980).  This  shallow  convergent  layer  (Fig.  4.1a)  in 
the  analyzed  case  is  much  greater  in  magnitude  than  the 
upper-level  outflow,  which  occurs  over  a  deeper  layer.  The 
absence  of  vertical  variations  in  the  LND  during  the  course 
of  this  storm  is  different  from  the  explosive  maritime 
cyclogenesis  case  studied  by  Calland  (1983),  in  which  the 
convergent  layer  was  found  to  deepen  as  the  storm  developed. 
Possible  explanations  for  this  are  the  lack  of  a  clearly 
defined  explosive  cyclogenesis  phase,  stronger  frictional 
effects  occur  over  land,  and  the  better  data  coverage  over 
land  which  provides  a  more  accurate  vertical  structure  of 
the  atmosphere. 

Much  larger  values  are  observed  for  the  inward  mass 
transport  in  the  analyzed  case  near  the  last  time  period 
(Fig.  4.1a).  This  suggests  that  some  mesoscale  feature  not 
captured  within  the  4°  radius  budget  volume  during  the 
earlier  time  periods  has  moved  into  the  budget  volume  at 
this  last  time  period.  An  oscillating  trend  in  the  forecast 
field  is  seen  with  maxima  at  12  GMT  28  March  and  06  GMT  29 
March  1984.  These  maxima  of  inward  horizontal  transport 
occur  about  6  h  prior  to  the  peak  model  precipitation 
periods  of  18  GMT  28  March  and  12  GMT  29  March  1984.  This 
strong  oscillation  pattern  does  not  occur  at  radii  greater 
than  4^  lat. 

The  area-averaged  kinematic  vertical  velocity  for  the 
analyses  (Fig.  4.2a)  indicates  an  absolute  maximum  near  700 
mb  for  the  00-12  GMT  29  March  1984  time  period.  The  fore¬ 
cast  area-averaged  omega  field  (Fig.  4.2b)  displays  an  abso¬ 
lute  maximum  at  03  GMT  29  March  1984  and  a  secondary  maximum 
from  09-21  GMT  28  March  1984  near  500  mb.  The  differences 
in  the  level  of  maximum  omega  are  sensitive  to  the  vertical 
profile  of  diabatic  heating  (Gyakum,  1983b).  Strong 
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diabatic  heating  in  the  lower  troposphere  is  believed  to 
induce  greater  vertical  velocities  and  thus  increase  the 
low-level  convergence  which  is  necessary  for  the  deepening 
of  the  cyclone.  In  the  forecast  case  (Fig.  4.2b),  there  is 
a  slight  bias  towards  more  intense  vertical  velocities.  The 
two  distinct  maxima  of  omega  evident  in  the  forecast  case 
will  have  important  implications  for  the  model-predicted 
precipitation  pattern  to  be  discussed  in  the  next  chapter. 

C.  VORTICITY  BUDGET 

In  the  vorticity  budget,  the  goal  is  to  determine  the 
relative  magnitudes  of  the  terms  in  the  vorticity  equation, 

6^a/6t  =  -1/A  $  Ca(Vn-  VQ)  dl  -d/dp  (*Si£) 

-  -  _  _  (4.1) 

-Ca  (7  *  v)  +  k  *  ( dV /dp  x  ?w)  +  Fr  +  R  , 


where  £a  is  the  absolute  vorticity,  A  is  the  area  of  a 
particular  surface  within  the  budget  volume,  Vn  is  the 


normal  wind  component,  VQ  is  the  cyclone's  normal  velocity 


component,  to  is  the  vertical  velocity,  Fr  denotes  friction, 
R  is  the  residual  and  the  overbar  denotes  an  area-average  of 
a  term.  The  term  on  the  left  side  of  (4.1)  is  the 
quasi-Lagrangian  tendency  of  vorticity.  On  the  right  side 
of  (4.1),  the  first  term  is  the  horizontal  transport  (flux) 
of  vorticity,  while  the  second  term  is  the  vertical  trans¬ 
port.  The  vorticity  divergence  term  is  the  third  term  on 
the  right  side  of  (4.1)  and  is  considered  as  a  source  or 
forcing  term  which  generates  positive  vorticity.  The  fourth 
term  on  the  right  side  is  the  tilting  term  and  the  fifth 
term  is  the  friction  or  dissipation  term.  The  residual  is 
the  last  term  on  the  right  side  of  (4. 1)  and  includes  the 
sources  and  sinks  not  explicity  resolved  by  the  other  terms 
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and  the  effects  of  computational  errors.  If  all  the  sources 
and  sinks  contributing  to  vorticity  changes  have  been  prop¬ 
erly  accounted  for  and  interpolation  errors  have  not  be 
introduced,  the  size  of  the  residual  should  be  relatively 
small.  Interpolation  and  calculation  errors  are  the  major 
source  of  uncertainty  in  budget  studies  based  on  numerical 
model  outputs.  For  the  analyses,  an  uncertainty  due  to 
incomplete  observations  is  also  present. 

The  lateral  transport  of  absolute  vorticity  can  be 
partitioned  into  mean  and  eddy  modes,  which  may  be  referred 
to  as  symmetric  and  asymmetric  components.  The  mean  mode 
represents  the  effects  due  to  mean  cyclone  convergence  and 
divergence,  whereas  the  eddy  mode  represents  the  horizontal 
transport  due  to  asymmetries  in  the  cyclonic  flow  which  lead 
to  correlations  between  the  wind  field  and  the  vorticity 
deviations.  Jet  streaks  and  short  wave  troughs  are  two 
mechanisms  which  contribute  to  the  eddy  mode  component. 

The  vertical  component  of  vorticity  is  computed  from  the 
horizontal  wind  field  by  using  a  finite  difference  form  of 
relative  vorticity  in  cylindrical  coordinates  and  adding  the 
appropriate  value  of  the  Coriolis  parameter  at  the  grid 
point.  The  NORAPS  forecast  and  analyzed  winds  are  greatly 
smoothed  by  passing  a  25-point  filter  eight  times  for 
display  purposes.  This  smoothing  can  also  lead  to  accumu¬ 
lated  errors  in  the  budget  calculations.  The  absolute 
vorticity  budget  results  for  the  NORAPS  forecast  and  anal¬ 
yses  are  presented  by  vertical/time  sections  for  radius  4° 
lat. 

D.  VORTICITY  BUDGET  RESULTS 

The  time  tendency  of  absolute  vorticity  is  evaluated 
using  a  forward  time  difference.  Since  this  finite  differ¬ 
ence  approximation  yields  a  time  tendency  at  the  mid-point 
of  the  two  time  periods,  the  remainder  of  the  terms  are 
averaged  over  two  time  periods  to  present  all  the  terms  in 
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the  budget  at  common  times.  For  the  forecast  case,  15  GMT 
28  March  refers  to  the  12-18  GMT  28  March  1984  time  period. 
The  analyses  times  refer  to  the  time  period  between  12  h 
synoptic  times  (00  and  12  GMT). 

The  time  tendencies  of  the  analyzed  vorticity  (  Fig. 
4.3a)  feature  a  large  positive  value  at  275  mb  for  the 
initial  period,  a  weak  tendency  pattern  at  the  middle 
period,  and  a  slightly  increasing  pattern  in  the  lower 
troposphere  for  the  last  time  period.  The  early  increase  at 
the  upper  levels  can  be  related  to  the  intensifying  jet 
streak  moving  into  the  western  part  of  the  budget  volume 
(Fig.  3.2).  The  forecast  time  tendency  of  absolute 
vorticity  (Fig.  4.3b)  reveals  the  same  general  pattern  as 
the  analyzed  case,  although  with  slightly  greater  magni¬ 
tudes.  The  negative  time  tendencies  in  the  middle  and  upper 
levels  during  the  18  GMT  28  March  time  frame  are  difficult 
to  explain  as  the  storm  circulation  intensifies.  Calland 
(1983)  found  a  similar  feature  in  a  North  Pacific  case  study 
and  suggested  this  decrease  was  due  to  a  subsynoptic  scale 
system  in  the  wind  field  during  that  period.  After  18  GMT 
28  March  1984,  positive  tendencies  occur  in  both  the  model 
and  analyzed  troposphere  as  the  increased  inward  horizontal 
mass  transport  (convergence)  spins  up  the  cyclone. 

The  lateral  transport  of  vorticity  (Fig.  4.4a)  in  the 
analyses  depicts  a  strong  inward  transport  of  vorticity 
below  775  mb  and  a  moderate  inward  transport  in  the  upper 
levels  of  the  troposphere.  A  relative  minimum  of  inward 
transport  is  located  at  600  mb.  The  transport  of  vorticity 
oscillates  in  the  vertical  between  inward  and  outward  values 
from  00-06  GMT  29  March  1984.  The  horizontal  transport 
(Fig.  4.4b)  for  the  forecast  case  verifies  well  in  the 
lower  troposphere  as  the  same  general  structure  is  depicted, 
although  somewhat  weaker.  However,  an  anomalous  feature  of 
outward  vorticity  transport  occurs  between  200  and  600  mb 
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from  09-21  GMT  28  March  1984  in.  the  forecast  case,  which  is 
not  verified  by  the  analyses.  The  reason  for  this  is  the 
large  divergence  contribution  which  overcompensates  the 
advective  contribution  during  this  time  interval.  Rather, 
an  inward  transport  of  vorticity  at  upper  levels  persists 
through  this  period  in  the  analyses. 

The  vertical  transport  of  vorticity  serves  to  vertically 
redistribute  the  vorticity  brought  into  the  budget  volume  by 
the  lateral  transport.  In  the  analyses  (Fig.  4.5a),  low- 
level  vorticity  is  exported  aloft,  resulting  in  weak 
vorticity  increases  in  the  middle  troposphere.  Above  350 
mb,  a  large  amount  of  vorticity  is  fluxed  vertically  into 
the  budget  volume  associated  with  the  strong  winds  at  the 
jet  stream  level  which  produce  large  vertical  gradients  of 
vorticity.  Another  possible  explanation  for  the  extremely 
high  values  in  the  upper  levels  of  the  model  is  that  a 
smaller  pressure  interval  (50  mb)  can  introduce  significant 
errors  in  conjunction  with  large  wind  speeds  when  the  finite 
difference  of  the  flux  is  performed.  The  forecast  case 
(Fig.  4.5b)  verifies  quite  well  with  a  slightly  higher  level 
(600  mo)  where  the  flux  changes  sign  from  upward  to  inward, 
while  slightly  weaker  features  are  predicted  above  600  mb. 

The  divergence  term  in  the  vorticity  budget  equation  is 
often  viewed  as  a  forcing  term  which  spins  up  the  vortex 
through  low-level  convergence.  A  relatively  deep  layer  of 
vorticity  convergence  (  to  350  mb)  is  seen  in  the  analyses 
(Fig.  4.6a)  until  00  GMT  29  March  1984.  This  layer  of 
vorticity  convergence  shrinks  by  300  mb  at  00  GMT  29  March 
1984.  This  may  indicate  a  separate,  distinct  smaller-scale 
and  shallower  low  center  not  previously  being  followed  has 
entered  the  budget  volume  and  results  in  the  dramatic 
increase  in  vorticity  convergence  below  700  mb.  By 
contrast,  the  forecast  case  starts  with  a  shallower  conver¬ 
gence  layer  which  deepens  throughout  the  forecast  period  to 
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450  mb.  The  distinct  maximum  in  vorticity  convergence  below 
700  mb  verifies  well.  A  forecast  maximum  in  upper-level 
divergence  occurs  at  21  GMT  28  March  1984,  which  is  9  h 
earlier  than  in  the  analyses. 

The  tilting  term  arises  from  the  vertical  velocity 
components  generated  when  horizontally  oriented  vorticity 
elements  are  tilted  to  the  vertical  by  a  horizontally 
varying  vertical  motion  field.  Thus,  the  vertical  shear 
tilts  the  horizontal  vorticity  elements  into  the  vertical 
position.  The  tilting  term  (Fig.  4.7a)  in  the  analyses  has 
a  weak  negative  contribution  throughout  the  depth  of  the 
atmosphere.  For  the  forecast  case  (Fig.  4.7b),  the  same 
general  pattern  of  weak  negative  values  verifies  in  the 
middle  troposphere  with  weak  positive  contributions  near  the 
top  and  bottom  of  the  model  atmosphere.  By  contrast, 
Calland  (1983)  found  the  tilting  term  to  be  a  source  of 
vorticity  in  his  case  study. 

Frictional  dissipiation  is  assumed  to  occur  only  in  the 
lowest  layer  of  the  model  and  is  parameterized  using  a 
stability  dependent  scheme  (Johnson  and  Downey,  1975).  The 
diagnosed  friction  (Fig.  4.8a)  plays  a  very  minor  role 
until  18  GMT  28  March  1984  when  the  circulation  has  become 
well  organized  and  leads  to  frictional  effects  that  oppose 
further  increases  in  vorticity.  This  frictional  dissipation 
is  strong  at  the  inner  radii  (less  than  6°  lat.  radius), 
which  suggests  that  friction  is  a  major  contributor  to  the 
vorticity  budget  only  near  the  central  region  of  high  winds. 
In  the  forecast  case,  the  frictional  dissipation  (Fig.  4.8b) 
verifies  well  with  slightly  greater  values  at  the  last  time 
period. 

The  residual  term  contains  accumulated  errors  resulting 
from  the  budget  calculations  and  any  sources  or  sinks  of 
vorticity  not  explicitly  computed.  A  positive  (negative) 
residual  in  the  vorticity  budget  indicates  an  apparent 


vorticity  source  ( sink)  in  which  the  observed  increases 
(decreases)  at  a  point  are  larger  than  is  estimated  by  the 
computed  terms.  Sources  can  be  due  to  such  things  as 
cumulus  friction  which  serves  to  transfer  low-level 
vorticity  aloft.  An  incorrect  representation  of  friction 
can  lead  to  a  sink  of  vorticity.  Residuals  for  forecast 
fields  should  be  nearly  zero  if  the  sources  and  sinks  of 
vorticity  have  been  properly  accounted  for  in  a  finite 
difference  representation.  Inaccuracies  in  vhe  vertical 
motion  and  horizontal  wind  analyses  contribute  to  the  phys¬ 
ical  errors  in  the  residual,  while  spatial  and  temporal 
finite  differencing  and  the  interpolation  to  the  budget 
volume  are  causes  for  computational  error. 

In  the  analyses,  the  residual  (Fig.  4.9a)  depicts  large 
negative  values  through  the  depth  of  the  atmosphere  with  the 
exception  of  a  positive  region  between  350  and  700  mb  after 
18  GMT  28  March  1984.  The  residual  (Fig.  4.9b)  in  the  fore¬ 
cast  is  generally  smaller  than  the  analyses.  The  large 
negative  residual  at  low  levels  indicates  the  NORAPS  model 
parameterizes  friction  much  stronger  than  the  stability- 
dependent  scheme  used  in  the  budget  formulation.  The  large 
vorticity  sink  in  the  analyses  (Fig.  4.9a)  in  the  upper 
levels  at  18  GMT  28  March  1984  is  hardly  evident  in  the 
forecast  case. 

The  terms  of  the  vorticity  budget  are  vertically  aver¬ 
aged  over  the  1000-500  mb  layer  to  obtain  a  clearer  under¬ 
standing  of  the  contribution  of  each  term  in  the  spin-up  of 
the  low-level  vortex  for  this  relatively  shallow  cyclone 
case.  These  results  are  shown  for  the  analyses  (Fig.  4.10) 
and  the  forecast  (Fig.  4.11).  The  friction  term  is  not 
shown  as  it  plays  a  relatively  minor  role.  The  leading  term 
is  the  vorticity  divergence  term  for  both  the  analyses  and 
forecast.  The  lateral  transport  of  vorticity  is  normally 
the  second  leading  term,  except  at  18  GMT  28  March  1984  in 


the  analyses  (Fig.  4.10)  where  the  residual  becomes  as  large 
as  the  leading  term.  The  time  tendency  of  vorticity  and  the 
residual  generally  rank  as  the  third  or  fourth  leading  term 
in  the  budget.  The  small  magnitude  of  the  residual  in 
comparison  to  the  other  terms  indicates  a  favorable  budget 
result.  The  vertical  transport  and  tilting  term  are  also 
generally  much  smaller  in  magnitude  than  the  leading  terms. 
Therefore,  the  vorticity  divergence  term  serves  as  the  major 
source  of  vorticity  to  the  budget  volume  through  strong 
low-level  convergence. 

E.  SUMMARY 

In  summary,  the  mass  budget  findings  reveal  a  shallow 
layer  of  low-level  convergence  beneath  a  deeper  layer  of 
upper-level  divergence.  The  level  of  maximum  vertical 
velocity  appears  fairly  constant  ( 700  mb)  during  the  course 
of  the  storm.  The  lateral  transport  is  the  leading  term  in 
the  vorticity  budget.  A  linkage  exists  between  the  mass  and 
vorticity  budgets  as  a  consequence  of  the  vorticity  diver¬ 
gence  term.  Low-level  convergence  acts  to  spin-up  the  low- 
level  vortex  as  viewed  from  the  time  tendency  term  in  the 
vorticity  budget  results.  This  increase  in  vorticity  is 
aided  by  a  favorable  upper-level  divergence  pattern,  which 
serves  to  export  mass  aloft  and  strengthens  the  secondary 
circulation  ( in-up-out)  in  the  vertical  plane.  The  eddy 
lateral  transport  provides  a  large  inward  transport  of 
vorticity  above  700  mb  to  offset  the  relatively  minor 
contribution  from  the  mean  or  divergent  mode. 

The  vertical  transport  term  serves  to  transport 
vorticity  that  has  been  "spun-up"  at  low  levels  into  the 
middle  troposphere.  The  tilting  term  acts  as  a  vorticity 
sink,  which  is  inconsistent  with  earlier  budget  findings 
(Calland,  1983).  Friction  only  plays  a  significant  role  at 
the  inner  core  as  it  offsets  the  low-level  vorticity  spin-up 
during  the  later  stages  of  the  storm.  The  diagnosed 


residual  from  the  vorticity  budget  is  observed  to  be  the 
balance  of  the  lateral  transport  and  time  tendency  terms. 

The  following  chapter  will  discuss  heat  and  moisture 
budget  results.  The  cooperation  between  the  mass  and 
vorticity  fields  is  an  essential  factor  in  organizing  the 
heating  and  moistening  in  the  storm  environment.  The 
heating  and  moistening  effects  can  alter  the  vertical  struc¬ 
ture  of  the  vertical  velocity  profile  and  produce  feedback 
effects  which  change  the  mass  and  vorticity  balance. 
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Fig.  4.1  Horizontal  mass  flux  for  radius  4  lat.  in  the  (A) 
Analysis  and  (B)  Forecast.  Negative/positive  values  denote 
mass  flux  into/out  of  budget  volume.  Contours  are  5  x 
10**- 11  g/s-100  mb. 
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Fig.  4.2  Area-averaged  vertical  velocities  in  the  (A) 
Analysis  and  (B)  Forecast  for  radius  4.  Negative  val 
indicate  upward  vertical  motion.  Contours  are  5  x  10* 
mb/s. 
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3_  Area-averaged  time  tendency  of  vorticity  for 
in  the  (A)  Analysis  and  (B)  Forecast.  Solid/dashed 
i vorticity  increases/decreases.  Contours 
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Fig.  4.10  Lower  troposDhere  (  1000-500  mb)  averacre  vorticitv 
budget  terms  for  analyses  at  radius  4:  Latera’l  transport 
(dashed),  time  tendency  (solid),  vorticity  divergence 
(dotted),  vertical  transport  (chain-dotted),  residual 
(  chain-asteri rk )  and  tilting  term  (chain-dashed).  Units  are 
1  x  10**-11  s*  * -2 . 


A.  GENERAL 


The  objective  of  performing  heat  and  moisture  budgets 
using  NORAPS  data  is  to  better  understand  the  physical 
processes  of  diabatic  heating  and  moistening  in  the  evolu¬ 
tion  of  explosive  maritime  cyclones.  Specifically,  we  want 
to  determine  the  leading  terms  of  each  equation  and  vertical 
distribution  of  heating  and  moistening  in  the  troposphere. 
The  diabatic  heating  and  moistening  rates  are  only  available 
as  a  residual  in  the  budgets  based  on  the  analysis  fields. 
However,  a  direct  comparison  can  be  made  between  the  rates 
diagnosed  from  the  forecast  fields  in  the  budget  formulation 
and  the  actual  rates  predicted  in  the  model.  In  particular, 
a  direct  comparison  can  be  made  of  the  level  and  phase  of 
the  maximum  heating  and  moistening  rates. 

The  moisture  budget  is  only  performed  for  the  forecast 
fields  since  a  moisture  analysis  is  not  performed  by  NORAPS. 
The  moisture  source  (residual)  at  model  levels  is  vertically 
integrated  for  comparison  with  the  area-averaged  6  hourly 
precipitation  amounts  from  the  model.  The  heat  and  moisture 
equations  used  for  the  budget  calculations  are  presented  at 
the  beginning  of  their  respective  sections.  Vertical/time 
sections  for  the  various  terms  in  each  equation  will  be 
presented  for  radius  4^  lat.  at  the  end  of  the  chapter.  The 
results  at  this  radius  represent  the  inner  core  of  the 
cyclone  where  heating  and  moistening  tend  to  be  concen¬ 
trated.  Results  are  available  out  to  radius  10^  lat.  but 
not  shown  due  to  space  limitations.  After  a  comparison 
between  the  heating/moistening  rates  diagnosed  from  the 
budget  calculations  and  the  corresponding  model-predicted 
rates,  the  column-integrated  results  are  given. 


B.  HEAT  BUDGET  RESULTS 
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The  area-averaged  flux  form  of  the  thermodynamic  equa¬ 
tion  is  presented, 

«T/St  =  -1/A  $  T(Vn-V0)  dl  (5.  D 

-d/dp  ( toT )  +  wa/cp  +  Q  , 

where  T  is  the  temperature,  A  is  the  area  of  a  surface 
within  the  budget  volume,  V  is  the  normal  wind  component,  VQ 
is  the  cyclone's  normal  velocity  component,  to  is  the 
vertical  velocity,  a  is  the  specific  volume,  cp  is  the 
specific  heat  at  constant  pressure  and  an  overbar  denotes  an 
area-average  of  a  term.  The  term  of  the  left  side  of  (5.1) 
is  the  quasi-Lagrangian  temperature  tendency.  On  the  right- 
hand  side  (RHS)  of  (5.1),  the  first  two  terms  are  the  hori¬ 
zontal  and  vertical  heat  fluxes.  The  third  term  on  the  RHS 
is  the  energy  conversion  term  and  the  final  term  is  the 
residual.  The  residual  is  a  measure  of  the  diabatic  heating 
plus  computational  errors  in  the  budget  calculation.  Using 
the  flux  form  of  the  equation  ( as  opposed  to  the  advective 
form)  reduces  the  computational  errors  and  attaches  greater 
physical  significance  to  the  diabatic  heating  term,  which  is 
not  directly  observed.  A  vertical  average  can  be  applied  to 
each  term  of  (5.1)  to  obtain  column-averaged  heating  rates. 
This  representation  serves  to  highlight  the  important  phys¬ 
ical  processes  pertaining  to  heating  or  cooling  during  the 
growth  of  the  storm.  These  results  are  presented  at  the  end 
of  this  section. 

The  vertical  time  sections  of  each  term  in  the  thermody¬ 
namic  equation  is  presented  in  the  advective  form  of  the 
equation  to  more  easily  interpret  the  physical  processes. 
The  energy  conversion  and  the  vertical  advection  terms  are 
shown  individually  and  then  combined  to  form  the  adiabatic 
cooling  term. 
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In  the  quasi-Lagrangian  temperature  tendency  term  ( Fig. 
5.1a),  the  analyzed  case  depicts  a  pattern  of  early-period 
cooling,  mid-period  warming,  followed  by  late-period  cooling 
below  600  mb.  This  pattern  can  be  interpreted  by  viewing 
the  cyclone  track  in  Fig.  3. la.  The  analyzed  low  center  in 
eastern  Texas  is  initially  in  an  extremely  warm  environment. 
As  the  storm  moves  northeast  into  a  colder  environment,  the 
cyclone  experiences  cooling.  The  cyclone  subsequently 
tracks  eastward  and  the  self-amplification  process  (building 
of  the  downstream  ridge)  causes  a  warming  at  18  GMT  28  March 
1984  through  the  entire  troposphere.  At  06  GMT  29  March 
1984,  low-level  cooling  appears  as  the  storm  is  nearing  its 
mature  stage  and  lower  temperatures  have  filled  most  of  the 
cyclone  volume  below  700  mb.  The  continued  heating  above 
700  mb  is  a  feature  that  will  be  discussed  following  the 
presentation  of  the  moisture  budget  results.  The  tempera¬ 
ture  tendency  (Fig.  5.1b)  in  the  forecast  case  shows  the 
same  general  pattern  with  the  cooling  at  the  last  period 
extending  to  a  higher  level  (350  mb).  By  comparison  with 
the  analyzed  case,  the  magnitude  of  this  term  is  slightly 
greater  (2°C/day)  overall  for  the  forecast  case  with  sharp 
differences  between  175  and  275  mb  at  09  GMT  29  March  1984. 
Magnitudes  of  12^C/day  are  not  considered  physically  real¬ 
istic  in  this  layer. 

As  the  horizontal  advection  in  the  quasi-Lagrangian 
formulation  has  the  advection  effects  due  to  cyclone 
velocity  removed,  this  term  isolates  the  purely  develop¬ 
mental  effects  of  horizontal  temperature  advection  as  the 
storm  intensifies.  This  term  is  computed  by  subtracting  the 
velocity  convergence/divergence  multiplied  by  temperature 
from  the  total  horizontal  temperature  flux  divergence.  The 


horizontal  temperature  advection  (Fig.  5.2a)  in  the  analyses 
depicts  unrealistically  large  magnitudes  (20^C/day)  near  225 
mb.  If  equal  amounts  of  cold  and  warm  advection  occupied 
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the  budget  volume  then  small  magnitudes  of  this  term  would 
be  realized.  Cold  advection  in  the  early  periods  is  evident 
below  450  mb  which  rapidly  shallows  at  12  GMT  28  March  1984. 
This  is  consistent  with  the  discussion  in  Chapter  3,  wherein 
stronger  cold  advection  is  found  in  an  area-averaged  sense 
early  in  the  storm,  while  the  warm  advection  pattern  is 
relatively  weak.  The  horizontal  advection  (Fig.  5.2b)  in 
the  forecast  case  verifies  well  with  the  exception  of  higher 
values  ( 30-35®C/day )  near  225  mb. 

A  disturbing  factor  is  the  unusually  large  values  of 
horizontal  temperature  advection  near  225  mb  in  both  the 
forecast  and  analyzed  fields.  In  general,  one  would  expect 
to  see  weak  warm  advection  aloft  within  the  budget  volume 
due  to  the  westward  tilt  of  a  baroclinic  system.  Since  the 
budget  volume  translates  with  the  surface  low  center,  warm 
advection  aloft  occurring  to  the  east  of  the  upper-level 
trough  position  would  be  a  displaced  to  the  west  of  the 
budget  volume.  Therefore,  relatively  weak  upper-level  hori¬ 
zontal  temperature  advection  would  be  occurring  within  the 
budget  volume.  Such  large  values  at  this  level  are  clearly 
in  error.  A  possible  explanation  for  the  forecast  case  is 
that  large  errors  can  be  introduced  in  the  interpolation 
from  sigma  to  pressure  surfaces.  However,  since  this 
problem  likewise  occurs  in  the  analyses,  a  faulty  interpola¬ 
tion  scheme  in  a  region  of  few  observations  can  account  for 
this  problem.  An  alternate  explanation  is  related  to  super¬ 
saturating  the  upper  levels  of  the  atmosphere  which  releases 
large  amounts  of  latent  heat,  as  will  be  discussed  in  more 
detail  in  the  next  section.  However,  this  would  only  be 
true  for  the  forecast  case. 

The  vertical  temperature  advection  (Fig.  5.3a)  in  the 
analyses  has  a  maximum  at  350  mb  with  a  40^C/day  rate  at  06 
GMT  29  March  1984.  This  maximum  level  in  the  forecast  (Fig. 
5.3b)  verifies  well,  but  a  much  higher  value  (65^C/day)  is 


predicted.  As  the  area-averaged  levels  of  maximum  omega 
(Fig.  3.2)  are  700  and  500  mb  for  the  analyses  and  forecast 
fields,  clearly  a  large  vertical  temperature  gradient  causes 
the  higher  maximum  level  of  vertical  temperature  advection. 
Large  vertical  gradients  in  temperature  at  350  mb,  where  one 
would  normally  expect  to  find  the  tropopause  for  such  a 
vigorous  extratropical  cyclone,  would  not  be  realistic. 

The  area-averaged  temperature  field  (Fig.  5.4)  indicates 
a  decrease  in  temperature  through  the  depth  of  the  atmos¬ 
phere  (to  100  mb),  without  showing  a  tropopause  in  both  the 
forecast  and  analyzed  cases.  Individual  grid  points  were 
examined  to  confirm  these  unrealistic  vertical  temperature 
structures.  These  unrealistic  vertical  temperature  profiles 
introduce  errors  into  the  column-averaged  heat  budgets  at 
upper  levels.  Vertical  temperature  profiles  in  the  data 
should  be  closely  checked  to  ensure  physical  consistency 
prior  to  undertaking  budget  computations. 

The  level  where  energy  conversion  (Fig.  5.5a)  is  a 
maximum  ( -55°C/day)  is  between  350  to  400  mb.  The  forecast 
case  (Fig.  5.5b)  depicts  a  much  different  profile  with  an 
absolute  maximum  of  -90°C/day  near  350  mb  at  03  GMT  29  March 
1984.  A  secondary  maximum  of  -45^C/day  is  found  at  09  GMT 
28  March  1984  at  the  same  level.  The  general  pattern  of  the 
vertical  profiles  for  both  the  analyzed  and  forecast  cases 
looks  very  similar  to  the  vertical  velocity  profiles  in  Fig. 
3.2,  except  the  maxima  are  shifted  upward  due  to  the 
weighting  factor  of  specific  volume  in  the  energy  conversion 
term. 

To  more  easily  facilitate  the  determination  of  the  resi¬ 
dual  by  viewing  the  contribution  of  individual  terms,  the 
energy  conversion  and  vertical  temperature  advection  are 
grouped  together.  These  two  terms  have  a  common  factor  of 
omega  and  are  of  opposite  signs.  The  resulting  term  is 
known  as  the  adiabatic  cooling  term,  written  as  to(  I -  f) 
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where  is  the  dry  adiabatic  lapse  rate  and  T  is  the  envi¬ 
ronmental  lapse  rate  or  vertical  temperature  gradient. 

-  r  is  always  constrained  to  be  greater  than  zero  in  the 
model  atmsophere  (due  to  the  dry  convective  adjustment 
scheme).  As  the  mean  vertical  motion  is  upward  in  the 
region  of  the  cyclone,  i.e,  to  is  less  than  zero,  this  term 
is  negative  and  opposes  heating.  This  term  can  also  be 
viewed  as  modulating  the  effects  of  latent  heating  during 
cyclone  development.  The  adiabatic  cooling  profile  (Fig. 
5.6a)  in  the  analyses  depict  the  strongest  cooling  tendency 
( -25°C/day)  at  850  mb  during  the  last  time  period.  In  the 
forecast  case  (Fig.  5.6b),  a  noticeable  difference  is  seen 
above  350  mb  with  strong  cooling  rates  of  25  to  35°C/day 
which  are  10  to  15®C/day  greater  than  the  analyses.  The 
25®C/day  cooling  rate  at  03  GMT  29  March  1984  is  300  mb 
above  the  maximum  in  the  analysis. 

The  surface  sensible  heat  flux  (Fig.  5.7)  is  archived 
for  the  forecast,  but  is  not  available  for  verification 
purposes  since  it  is  not  directly  observable.  A  definite 
diurnal  trend  is  seen  in  the  surface  sensible  heat  flux 
pattern  because  most  of  the  storm  volume  is  over  land.  This 
term  contributes  to  the  warming  of  the  column  below  775  mb. 
In  a  column  and  area-averaged  sense,  its  overall  contribu¬ 
tion  is  very  small.  A  maximum  heating  rate  of  10®C/day 
occurs  from  15-21  GMT  28  March  1984  (daylight  hours)  below 
850  mb.  However,  its  effect  is  significant  over  localized 
regions  since  this  upward  flux  of  heat  can  destabilize  lapse 
rates  and  initiate  convection  and  be  an  important  contrib¬ 
utor  during  a  particular  phase  of  storm  development. 

The  residual  term  is  computed  as  the  balance  of  the 
other  calculated  terms,  and  thus  includes  diabatic  heating 
and  computational  errors.  Since  such  physical  processes  as 
latent  heating  and  radiational  cooling  cannot  be  directly 
measured,  they  must  be  inferred  as  a  residual  in  the 


computation.  The  radiational  cooling  term  is  generally 
acknowledged  to  be  small  in  relation  to  the  latent  heating. 
The  residual  in  the  budget  based  on  analyses  (Fig.  5.8a) 
depicts  the  level  of  maximum  heating  at  550  mb.  The  abso¬ 
lute  maximum  in  heating  (  20^C/day)  occurs  at  06  GMT  29  March 
1984.  This  heating  rate  is  in  close  agreement  with  the 
results  of  Liou  and  Elsberry  (1985),  who  used  European 
Center  for  Medium  Range  Weather  Forecasting  (ECMWF)  analyses 
for  a  western  North  Pacific  Ocean  case  study  and  found  a 
25°C/day  heating  rate  at  600  mb.  By  comparison  with  the 
residual  in  the  forecast  case  (Fig.  5.8b),  the  level  of 
maximum  heating  verifies  well  (500-550  mb)  and  a  secondary 
maximum  is  predicted  between  225  and  350  mb  which  is  not 
realistic.  Large-scale  precipitation  and  middle-level 
convection  occurring  within  the  extratropical  cyclone  would 
produce  a  lower  level  of  maximum  heating  as  opposed  to  a 
tropical  cyclone  where  deep  convection  would  cause  greater 
diabatic  heating  aloft.  Two  separate  maxima  in  heating 
occur  at  09  GMT  28  March  ( 15°C/day)  and  03  GMT  29  March  1984 
(20°C/day).  Another  difference  from  the  analyses  fields  is 
the  cooling  occurring  below  850  mb  at  03  GMT  29  March  1984, 
which  is  probably  due  to  evaporational  cooling. 

The  vertical  profile  of  the  heating  rates  (Fig.  5.9) 
predicted  by  the  model  is  available  for  comparison  with  the 
budget  diagnosed  residual  from  the  forecast  fields.  This 
heating  rate  is  the  average  of  two  instantaneous  values  each 
6  h  during  the  integration  and  is  not  the  accumulated  6  h 
total  which  would  be  more  comparable  to  the  values  derived 
diagnostically.  The  instantaneous  values  ( ^C/time  step) 
have  been  converted  to  a  daily  rate  (^C/day).  These  heating 
rates  include  net  condensational  heating  plus  short  and 
long-wave  radiational  effects.  The  surface  sensible  heat 
flux  is  an  output  field  of  the  NORAPS  model  but  is  not 
included  in  the  model  heating  rate  since  it  is  not 


distributed  over  a  constant  depth  in  the  model  prediction  as 
is  assumed  in  the  diagnostic  model  (lowest  200  mb).  The 
level  of  maximum  heating  is  450  mb  (Fig.  5.9a)  which  is 
higher  than  the  level  of  550  mb  diagnosed  in  the  analyses 
and  forecast  from  the  budget  residual  calculation.  Two 
distinct  maxima  appear  in  the  model-predicted  vertical 
heating  profile:  a  30^C/day  center  at  15  GMT  28  March  and  a 
35®C/day  center  at  03  GMT  29  March  1984.  Although  the  pres¬ 
ence  of  a  double  maximum  pattern  agrees  with  the  diagnosed 
residual  from  the  analyzed  fields,  the  magnitudes  of  the 
model-predicted  heating  rates  are  about  100  and  50%  greater, 
respectively.  Values  at  individual  grid  points  are  on  the 
order  of  200-400°C/day  in  localized  regions  of  the  storm. 

The  model-predicted  rates  are  also  approximately  double 
those  diagnosed  for  the  forecast  fields.  This  difference 
can  be  easily  shown  by  subtracting  the  model  heating  rates 
from  those  diagnosed  in  the  budget  for  the  forecast  fields 
(Fig.  5.9b).  This  figure  gives  us  an  estimate  of  calcula¬ 
tion  errors  in  the  budget.  Ideally,  we  would  desire  to  have 
very  small  differences  to  verify  the  accuracy  of  the  budget. 
The  consistent  profile  of  errors  reveals  an  overestimation 
of  heating  by  the  model  from  275  to  775  mb  and  excessive 
cooling  predicted  by  the  model  in  the  lowest  layers.  The 
errors  above  275  mb  are  a  result  of  the  spuriously  large 
values  of  horizontal  temperature  advection  in  the  budget 
calculation. 

The  column  and  area-averaged  heat  budget  results  are 
presented  for  the  analyses  (Fig.  5.10)  and  forecast  (Fig. 
5.11).  The  leading  term  (in  a  cooling  sense)  is  the  adia¬ 
batic  cooling  term,  which  is  strongly  modulated  by  the 
vertical  velocity  field.  The  diabatic  heating  term  ( resi¬ 
dual  plus  errors)  is  the  next  largest  term  and  is  on  the 
order  of  10  to  12^C/day  during  the  last  half  of  the  storm. 
The  primary  role  of  diabatic  heating  is  to  offset  adiabatic 


V 


w 


cooling  and  promote  heating  within  the  budget  volume  which 
can  contribute  to  storm  intensification.  The  horizontal 
temperature  advection  is  smaller  than  the  residual  during 
the  early  stage  of  the  storm  but  becomes  comparable  to  the 
residual  during  the  last  half  of  the  storm.  The 
quasi -Lagrangi an  temperature  tendency  contributes  only 
weakly  to  heating  or  cooling  the  budget  volume  during  the 
course  of  the  storm. 

The  near-zero  heating  rates  in  the  model  at  the  begin¬ 
ning  of  the  forecast  period  are  not  realistic  because  of  the 
extensive  convection  observed  along  the  Gulf  Coast.  The 
difference  between  the  area  and  column-averaged  budget  resi¬ 
dual  and  model-predicted  heating  rates  at  radius  4°  lat. 
(Fig.  5.12)  is  only  1  to  2°C/day  from  03  GMT  28  March  to  03 
GMT  29  March  1984.  The  quasi-Lagrangian  temperature 
tendency  can  introduce  errors  due  to  the  rapid  deceleration 
of  the  storm  during  the  last  12  h  of  the  model  storm. 
Although  not  previously  discussed,  the  column  and  area- 
averaged  residual  for  radius  6°  lat.  is  8.  5°C/day  during  the 
intensification  of  the  storm,  which  agrees  well  with  Bosse's 
( 1984)  results  for  the  same  radius  during  an  explosive  stage 
of  cyclogenesis  in  the  western  North  Pacific  Ocean.  Thus, 
v,e  are  provided  with  additional  evidence  that  this  storm  has 
many  similarities  to  an  explosive  maritime  cyclogenesis 
case. 

C.  MOISTURE  BUDGET 

An  important  topic  to  be  discussed  before  presenting 
moisture  budget  results  is  how  the  initial  moisture  field  is 
specified  in  the  NORAPS  model.  The  moisture  analysis  is 
obtained  by  using  the  12  h  NORAPS  forecast  specific  humidity 
field.  The  NORAPS  forecast  fields  of  moisture  content  are 
output  in  terms  of  vapor  pressure.  The  vapor  pressure  is 
internally  converted  into  specific  humidity  within  the 
budget  program. 
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As  with  the  other  budgets,  various  sources  of  errors  can 
arise.  The  interpolation  from  the  NORAPS  grid  to  the  budget 
volume  is  with  a  bi-cubic  spline  interpolation  scheme.  For 
very  small  values  of  vapor  pressure,  particularly  above  500 
mb,  small  negative  values  of  vapor  pressure  can  be  intro¬ 
duced.  Another  problem  can  arise  in  the  vertical  moisture 
flux  calculation.  A  harmonic-mean  scheme  is  used  in  NORAPS 
to  interpolate  the  moisture  field  (  at  odd  levels  in  the 
model)  to  the  even  levels  where  vertical  velocity  is  speci¬ 
fied.  This  scheme  has  the  effect  of  biasing  the  smaller 
values  downward,  which  can  cause  erroneously  large  moisture 
divergence/convergence  in  the  vertical.  Thus,  the  vertical 
flux  convergence  of  moisture  that  is  computed  in  the  budget 
model  by  using  a  linear  vertical  interpolation  scheme  can 
also  introduce  errors. 

A  leapfrog  (centered  differencing)  scheme  is  currently 
used  for  the  moisture  prediction  equation.  The  original 
version  of  NORAPS  used  a  forward-in-time  and  space  (upwind 
or  one-sided  difference)  scheme,  which  has  larger  truncation 
errors  and  tends  to  damp  the  smaller-scale  features  in  time. 
The  change  to  this  new  scheme  resulted  in  some  complica¬ 
tions.  One  problem  was  the  omission  of  a  time  filter  to 
smooth  the  large  computational  modes  that  cause  oscillations 
in  a  centered-dif f erencing  scheme.  The  resulting  problem 
for  the  budget  is  that  large  errors  can  occur  in  the 
quasi-Lagrangian  moisture  tendency  (either  too  small  or  too 
large ) . 

The  terms  of  the  area-averaged  moisture  budget  equation 
are  given  in  (5.2), 

6q/St  =  -  1/A  $  q(vn-vG)  dl 
-d/dp  (o>q)  +  E  -  P 


(5.2) 


where  q  is  the  specific  humidity,  A  is  the  area  of  a  partic¬ 
ular  surface  within  the  budget  volume,  Vn  is  the  normal  wind 
component,  VQ  is  the  cyclone's  normal  velocity  component,  to 
is  the  vertical  velocity,  E  is  evaporation  ,  P  is  precipita¬ 
tion  and  an  overbar  denotes  the  area-average  of  a  term.  The 
quasi-Lagrangian  moisture  tendency  is  the  term  on  the  left 
side,  whereas  the  horizontal  and  vertical  moisture  fluxes 
are  the  first  two  terms  on  the  right  side  of  (5.2),  respec¬ 
tively.  The  last  term  on  the  right  side  is  the  residual  or 
moisture  sink  term.  This  term  includes  the  combined  effects 
of  condensation,  evaporation  and  moisture  fluxes  from  the 
surface.  As  these  terms  are  not  directly  observable,  they 
must  be  computed  as  a  residual.  The  convention  for  this 
residual  is  negative  (positive)  if  precipitation  (evapora¬ 
tion)  is  occurring.  The  quasi-Lagrangian  moisture  tendency 
differs  from  the  moisture  sink  term  in  that  specific 
humidity  is  a  measure  of  the  amount  of  water  vapor  (relative 
humidity)  in  the  air,  whereas  the  moisture  sink  is  a  measure 
of  the  liquid  water  being  condensed  within  the  budget  volume 
(assuming  precipitation  exceeds  evaporation).  When  conden¬ 
sation  occurs  within  a  region,  the  phase  change  results  in  a 
decrease  in  the  amount  of  water  vapor.  In  practice,  while 
precipitation  is  measured  at  regular  intervals  (over  land) 
by  the  standard  observing  network,  surface  moisture  fluxes 
and  evaporation  are  not  directly  measured  but  must  be 
inferred  from  approximations. 

The  primary  objective  in  performing  the  moisture  budget 
is  to  determine  the  agreement  with  the  heat  budget  with 
respect  to  the  level  and  phase  of  maximum  heating  and  mois¬ 
tening.  As  in  the  other  budgets,  the  moisture  budget  for 
radius  4^  lat.  will  be  discussed.  Since  a  moisture  anal¬ 
ysis  is  not  performed  by  NORAPS,  a  comparison  of  the  diag¬ 
nosed  budget  results  in  the  forecast  case  cannot  be  made 
with  the  analyzed  case.  A  comparison  is  made  between  the 
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moisture  source  term  diagnosed  from  the  budget  and  the  mois¬ 
tening  rate  computed  directly  in  the  model.  The  moistening 
rate,  as  in  the  model-predicted  heating  rate,  is  the  average 
of  two  instantaneous  values  each  6  h  during  the  integration. 

D.  MOISTURE  BUDGET  RESULTS 

The  quasi-Lagrangian  moisture  tendency  (Fig.  5.13a)  in 
the  forecast  reveals  an  oscillating  trend  in  phase  with  the 
precipitation  cycle  (Fig.  5.14).  Relatively  strong  mois¬ 
tening  occurs  prior  to  09  GMT  28  March  1984,  and  is  followed 
by  a  weak  period  of  moisture  decrease  within  the  budget 
volume  from  09-18  GMT  28  March  1984.  This  period  of  mois¬ 
ture  decrease  correlates  with  the  precipitation  "burst" 
during  the  12-18  GMT  28  March  1984.  A  strong  increase  in 
the  moisture  tendency  centered  at  775  mb  occurs  from  18  GMT 

28  March  to  06  GMT  29  March  1984  prior  to  the  precipitation 
"burst"  between  06  and  12  GMT  29  March  1984  (Fig.  5.14).  A 
decrease  in  the  moisture  tendency  occurs  from  06  to  09  GMT 

29  March  1984  which  correlates  with  the  moisture  loss  from 
the  atmosphere  during  the  heavy  precipitation  period. 

The  horizontal  moisture  flux  (Fig.  5.15a)  derived  from 
the  forecast  fields  has  a  layer  of  moisture  convergence 
(  surface  to  500  mb)  beneath  a  shallower  and  much  weaker 
layer  of  moisture  divergence  from  275  to  500  mb.  Fuelberg 
et  al.  (1986)  observed  a  much  deeper  layer  (up  to  350  mb) 
of  moisture  convergence  in  a  severe  thunderstorm  environ¬ 
ment.  They  used  a  dense  observation  rawindsonde  network  in 
a  stationary  budget  study  over  a  relatively  short  time 
period  with  3  h  observations.  The  horizontal  moisture  flux 
is  composed  of  an  advective  and  convergence/divergence  term. 
The  convergence/divergence  term  accounts  for  75%  of  the 
total  moisture  flux  in  this  study,  which  is  similar  to  the 
results  of  Fuelberg  et  al.  (  1986).  The  reason  for  the 
weaker  moisture  convergence  aloft  in  this  study  is  attrib¬ 
uted  to  the  shallower  horizontal  convergence  layer  as  shown 
in  the  mass  budget  for  this  case. 
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The  vertical  moisture  flux  term  represents  the  redistri¬ 
bution  of  moisture  in  the  vertical.  Horizontal  moisture 
flux  convergence  at  the  lower  levels  must  experience  a 
corresponding  vertical  flux  divergence  from  continuity 
considerations.  The  vertical  moisture  flux  (Fig.  5.15b) 
based  on  the  forecast  fields  shows  the  pattern  of  low-level 
moisture  divergence  below  775  mb  and  moisture  convergence 
from  775  to  275  mb.  However,  an  alternating  and  unrealistic 
pattern  is  noted  in  the  vertical  profile  of  the  vertical 
moisture  flux  after  21  GMT  28  March  1984.  The  individual 
grid  points  were  examined  in  the  neighborhood  of  the  storm 
to  determine  which  variable  contributed  to  this  unrealistic 
profile.  The  400  mb  specific  humidity  field  at  06  GMT  29 
March  1984  (Table  1)  is  displayed  for  the  region  surrounding 
the  storm.  Unrealistic  values  in  excess  of  10  g/kg  are  seen 
in  the  fourth  and  fifth  rows  of  this  field.  Small  negative 
values  are  also  noted  in  this  field  which  may  be  due  to  the 
bi-cubic  horizontal  interpolation  with  such  large  moisture 
gradients.  An  oscillatory  trend  in  the  vertical  specific 
humidity  profile  (not  shown)  is  evident  between  850  and  400 
mb  at  06  GMT  29  March  1984,  which  is  directly  responsible 
for  the  vertical  flux  profile  (Fig.  5.15b). 

The  question  naturally  arises  as  to  the  cause  of  such 
unrealistic  specific  humidity  values.  Checks  (by  Dr.  C. -S. 
Liou)  of  the  NORAPS  program  revealed  two  sources  of  error. 
First,  an  incorrect  specification  of  the  moistening  rate 
occurred  during  the  switch  to  a  centered-dif ferencing 
scheme.  The  heating  rate  was  likewise  in  error  since  only 
one-half  the  appropriate  rate  was  used  in  the  heat 
prediction  equation.  The  moistening  rate  appropriate  for 
the  one-sided  differencing  scheme  was  used  in  the  centered- 
differencing  scheme,  and  therefore  was  too  low  by  a  factor 
of  two.  Since  the  (negative)  moistening  rate  is  only  one- 
half  of  the  true  precipitation  rate,  moisture  is 


continuously  carried  over  into  each  successive  time  step, 
which  leads  to  a  supersaturation.  Second,  the  omission  of  a 
time  filter  in  the  moisture  prediction  equation  using  the 
leapfrog  scheme  could  result  in  erroneous  specific  humidity 
values  if  left  unsmoothed. 

The  "over-moistening"  has  a  direct  effect  on  the  precip¬ 
itation  predicted  by  the  model.  The  model  precipitation 
field  (without  area-averaging)  has  a  maximum  value  exceeding 
10  cm/6  h  at  12  GMT  29  March  1984  (Fig.  5.14).  Maxima  in 
precipitation  during  the  36  h  forecast  period  occur  at  18 
GMT  28  March  and  12  GMT  29  March  1984.  A  nearly  linear 
increase  in  precipitation  for  this  storm  would  be  more 
consistent  with  available  observations  (over  land).  A 
review  of  NWS  daily  weather  summary  sheets  for  28  and  29 
March  1984  showed  the  largest  observed  rainfall  rates  (north 
of  the  warm  front)  to  be  on  the  order  of  5  cm/day. 
Area-averaged  precipitation  for  radius  4°  lat.  from  the 
NORAPS  model  ranges  between  2  and  4  cm/day  until  the  last 
time  period  when  the  rate  increases  to  6  cm/day. 

The  residual  in  the  moisture  budget  includes  the  conden¬ 
sation,  evaporation  and  the  surface  moisture  fluxes.  A 
negative  value  indicates  more  condensation  would  be  occur¬ 
ring  than  evaporation,  while  a  positive  value  indicates  more 
evaporation  would  be  occurring  than  condensation.  Since  no 
liquid  water  is  retained  in  the  model  atmosphere,  a  positive 
residual  at  upper  levels  indicates  an  unrealistic  result. 
The  vertical  cross  section  of  the  residual  (Fig.  5.16a)  in 
the  budget  based  on  forecast  fields  reveals  a  large  moisture 
sink  in  the  middle  troposphere  throughout  the  entire  period 
due  to  condensation  exceeding  evaporation.  The  maximum 
levels  of  moistening  occur  at  600  and  350  mb.  An  exception 
is  an  anomalous  positive  region  near  400  mb  at  03  GMT  29 
March  1984  that  is  associated  with  the  spurious  vertical 
flux  profile  discussed  above.  The  small  residuals  above  275 


mb  are  an  artifact  of  the  negligible  amount  of  moisture  in 
this  region.  An  absolute  maximum  in  moistening  occurs  at  09 
GMT  29  March  1984,  which  is  slightly  contaminted  by  the 
problem  discussed  with  the  profile  of  the  vertical  moisture 
flux.  This  strong  moisture  sink  (precipitation)  agrees 
favorably  with  the  large  precipitation  rates  predicted  by 
the  model  at  this  time.  Another  maximum  in  moistening 
occurs  at  15  GMT  28  March  1984  which  correlates  well  with 
the  precipitation  maximum  at  18  GMT  28  March  1984  ( Fig. 
5.  14). 

The  archived  moistening  rate  in  the  forecast  model  (Fig. 
5. 16b)  has  a  different  structure  with  much  greater  rates 
than  diagnosed  from  the  moisture  budget.  Distinctive  maxima 
are  located  near  450  mb  at  15  GMT  28  March  and  at  03  GMT  29 
March  1984.  These  moistening  rates  necessarily  correspond 
in  level  and  phase  to  the  heating  rates  (Fig.  5.9).  These 
archived  rates  are  approximately  double  those  diagnosed  in 
the  budget  calculation.  This  is  a  result  of  the  incorrect 
specification  of  the  moistening  rate  when  switching  to  a 
centered-difference  scheme,  which  introduces  a  factor  of  two 
difference.  The  column-integrated  moisture  budget  results 
for  the  diagnosed  residual  (moisture  source),  the  model- 
predicted  moistening  rate  and  the  area-averaged  precipita¬ 
tion  rates  are  presented  in  Fig.  5.  17.  As  the  moistening 
rate  predicted  by  the  model  contains  the  accumulated  effects 
of  condensation,  evaporation  and  surface  moisture  fluxes  of 
the  time  average  of  two  instantaneous  values  each  6  h  apart, 
it  differs  from  the  model  precipitation  which  is  an  accumu¬ 
lated  6  h  total.  In  fact,  the  two  would  be  exactly  equal  if 
an  accumulated  6-h  archived  moistening  rate  was  used  and 
there  were  no  computation  errors  in  the  budget.  However, 
the  diagnosed  moisture  source  from  the  budget  is  systemati¬ 
cally  less  than  the  model-predicted  precipitation  and  mois¬ 
tening  rates.  Since  this  error  is  systematic,  this  provides 
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further  evidence  of  a  model  problem  if  the  budget  calcula¬ 
tions  have  been  performed  properly.  This  difference  between 
the  curves  (Fig.  5.17)  can  be  traced  to  the  application  of 
only  one-half  the  moistening  rate  in  the  centered- 
differencing  scheme  for  the  moisture  prediction  equation. 
This  error  in  the  model  leads  to  an  overprediction  of  the 
moistening  and  precipitation  rates. 

E.  SUMMARY 

The  moisture  budget  results  show  the  pattern  of  horizontal 
moisture  flux  convergence  in  the  surface  layers  which  is 
modulated  by  the  vertical  velocity  field.  The  moisture  is 
then  transported  into  the  upper  troposphere  through  the 
vertical  moisture  flux,  cools  moist  adiabatically  until  it 
reaches  saturation  and  then  condenses  and  falls  as 
precipitation.  The  phase  change  from  water  vapor  to  liquid 
releases  heat  at  a  certain  level  in  the  troposphere  which  is 
an  essential  ingredient  to  the  energetics  of  the  cyclone. 

The  major  factor  for  many  of  the  errors  in  the  vertical 
distribution  of  heating  and  moistening  is  believed  to  be 
linked  to  the  incorrect  specification  of  the  heating  and 
moistening  rate  which  leads  to  an  "over-moistening"  of  the 
model  atmosphere.  The  harmonic-mean  scheme  used  in  calcu¬ 
lating  the  vertical  moisture  fluxes  for  the  moisture 
prediction  may  also  account  for  some  of  the  variance  in  che 
vertical  moisture  flux  during  the  last  two  time  periods 
(Fig.  5.15b). 

The  moisture  budget  results  reveal  a  discrepancy  between 
the  diagnosed  and  the  archived  precipitation  results.  The 
precipitation  fields,  the  vertical  moisture  flux  and  the 
model-predicted  moistening  rates  provide  evidence  of  a 
systematic  error  in  the  NORAPS  model.  A  distinct  oscil¬ 
lating  pattern  is  found  in  the  model  precipitation  field. 
This  oscillatory  trend  is  likely  linked  to  the  area-averaged 
vertical  motion  field  in  the  model  forecast  (Fig.  3.2b) 
which  also  exhibits  a  double  maximum  in  the  field. 
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Fig.  5.10  Area  and  column-averaged  heat  budget  terms  at 
raaius  4  based  on  analyzed  fields.  Solid  line  depicts 
temoerature  tendency  term,  dotted  line  is  adiabatic  cooling, 
dashed  line  is  horizontal  temperature  advection  and  tne 
chain-dotted  line  is  the  diabatic  heating.  Units  of  heating 
are  C/'day  and  the  time  is  given  in  day  (first  two  digits) 
and  hours  (in  Greenwich,  second  two  digits). 
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NORAPS  SPECIFIC  HUMIDITY  (G/KG)  VALUES  AT  400  MB  FOR 
06  GMT  29  MARCH  1984  (L  MARKS  POSITION  OF  THE  LOW) 
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Fig.  5.16  (A)  Residual  for  moisture  budqet  and  (B)  Model- 
predicted  moistening  rate  for  radius  4.  Contour  interval  is 
10  x  10**-4  g/g/day. 
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Fig.  5.17  Area-averaged,  column- intearated  moisture  budget 
results  for  radius  4.  Solid  line  represents  the  budget 
residual,  dashed  line  is  the  model-predicted  moistening  rate 
and  the  dotted  line  is  the  model  precipitation.  Units  are 
cm/day. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  physical  processes  in  the  development  of  the  extra- 
tropical  cyclone  occur  over  a  broad  spectrum  of  space  and 
time  scales.  This  study  has  concentrated  on  the  dynamic  and 
thermodynamic  forcing  on  the  cyclone  scale  as  it  relates  to 
storm  intensification.  The  low-level  mass  convergence 
serves  to  spin  up  the  low-level  vortex  as  viewed  in  the 
vorticity  budget  results.  Increased  upward  vertical  motion 
results  in  stronger  convergence,  which  serves  to  increase 
the  horizontal  moisture  transport  at  the  lower  levels.  This 
moisture  is  transported  aloft  as  seen  in  the  vertical  mois¬ 
ture  flux.  The  net  increase  in  moisture  within  the  column 
plus  the  cooling  process  (moist  adiabatic)  leads  to  precipi¬ 
tation,  which  greatly  exceeds  the  quasi-Lagrangian  moisture 
tendency.  The  release  of  latent  heat  contributes  to  a  posi¬ 
tive  quasi-Lagrangian  temperature  tendency  aloft  (Fig. 
5.1a),  while  cooling  occurs  in  the  lower  troposphere  in 
response  to  cold  advection,  evaporation  and  forced  lifting 
of  sub-cloud  air  by  diabatic  processes  in  the  cloud. 

In  particular,  special  emphasis  has  been  placed  on 
heating  and  moistening  which  are  essential  components  in  the 
cyclone  evolution.  Results  indicate  the  heating  and  mois¬ 
tening  rates  predicted  by  the  model  are  50  to  100%  greater 
than  the  rates  diagnosed  from  budget  residual  calculations. 
This  difference  is  largely  systematic  and  can  be  attributed 
to  the  incorrect  specification  of  the  moistening  and  heating 
rates  within  the  NORAPS  model.  Consequently,  the  heat  and 
moisture  budgets  suffer  in  terms  of  yielding  physically 
inconsistent  results  at  various  times  and  levels.  Revised 
heat  and  moisture  budget  results  (based  on  corrections  in 
predicting  the  heating  and  moistening  rates)  are  presented 
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in  Appendix  C  and  constitute  an  important  section  of  this 
study. 

The  release  of  latent  heat  caused  by  strong  inward  mois¬ 
ture  transport  in  the  lower  troposphere  into  the  storm's 
central  core  is  believed  to  play  a  critical  role  in  early 
intensification  of  the  extratropical  cyclone.  The  phase  and 
shallow  level  of  diabatic  heating  are  key  ingredients  in 
modulating  the  sea- level  pressure  tendency.  Previous  inves¬ 
tigations  by  Danard  (1966),  Gyakum  (1983b)  and  Anthes  et  al. 
(1983)  stress  that  diabatic  heating  only  contributes  to 
storm  intensification  during  the  later  stages.  This  view  is 
generally  accepted  since  a  sizeable  precipitation  field  is 
rarely  evident  in  the  early  stages  of  the  storm.  However, 
strong  conditional  instability  is  often  observed  upstream  of 
surface  troughs  during  the  initial  stages  of  baroclinic 
instability  over  ocean  regions.  As  a  consequence  of  the 
relatively  warm  ocean  surface  underlying  the  extremely  cold 
air  above  the  planetary  boundary  layer,  near-neutral  lapse 
rates  with  strong  low-level  moisture  convergence  cause  the 
atmosphere  to  respond  with  cumulus  convection.  The  latent 
heat  released  by  the  convection  along  with  cumulus- induced 
subsidence  play  a  major  role  in  strongly  heating  the  lower 
troposphere.  Chang  et  al.  (1982)  suggest  that  the  release 
of  latent  heat  aids  in  the  conversion  of  potential  to 
kinetic  energy  by  acting  as  a  catalyst.  Thus,  the  latent 
heating  in  a  favored  region  of  the  developing  storm  ( north¬ 
west  quadrant)  could  accelerate  the  the  baroclinic  insta¬ 
bility  process. 

There  is  a  definite  need  to  understand  how  diabatic 
heating  alters  the  mass  field  around  the  storm,  creating 
strong  dynamic  mesoscale  features.  Secondary  ageostrophic 
circulations,  forced  by  diabatic  heating,  have  been 
suggested  as  one  mechanism  where  the  mass  field  is  forced  to 
respond.  If  diabatic  heating  is  a  major  contributor  to  the 
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forcing  of  mesoscale  features,  which  are  believed  to  trigger 
explosive  cyclogenesis  events,  a  better  definition  of  these 
events  will  need  to  be  obtained  through  a  higher  resolution 
observational  network. 

B.  RECOMMENDATIONS 

The  use  of  the  higher  spatial  (horizontal  and  vertical) 
resolution  NORAPS  model  is  important  in  modeling  the  evolu¬ 
tion  of  mesoscale  features  which  affect  the  cyclone  scale. 
These  budget  results,  although  seemingly  contaminated  by 
several  sources  of  errors,  do  provide  a  baseline  for  future 
studies  using  the  NORAPS  model. 

The  following  recommendations  are  made  for  future 
studies: 

•  Resolve  problems  with  the  NORAPS  heating  and  moistening 
rates  for  a  more  accurate  comparison  of  budget  resi¬ 
duals  with  the  model-predicted  rates,* 

•  Initiate  a  moisture  analysis  for  NORAPS  to  better 
predict  precipitation  and  other  moistening  processes 
within  the  storm  environment. 

•  Use  higher  resolution  (spatial  and  temporal)  analyses 
to  validate  the  thermal  and  moisture  structure  of 
explosive  maritime  cyclogenesis  cases. 


NAVY  OPERATIONAL  REGIONAL  ATMOSPHERIC  PREDICTION  SYSTEM 


( NO RAPS) 

1.  MODEL  CHARACTERISTICS 

|  NORAPS,  which  includes  an  analysis  and  a  regional  fore- 

> 

j  cast  model,  produces  high  spatial  resolution  (typically  80 

km  in  the  western  Atlantic  version),  short  term  (36-48  hr) 
numerical  forecasts  over  a  limited  domain.  NORAPS  was 
developed  by  Dr.  Rich  Hodur  of  the  Naval  Environmental 
Prediction  Research  Facility  (NEPRF),  who  kindly  provided 
the  fields  used  in  this  research.  The  principle  advantage 
of  using  NORAPS  as  opposed  to  a  global  or  hemispheric  model 
is  the  small  spatial  scale  features  that  are  resolved.  This 
model  has  the  additional  asset  of  flexibility  as  the  grid  is 
globally  relocatable,  and  the  user  may  specify  the  dimen¬ 
sions  and  horizontal/vertical  resolution.  An  additional 
feature  is  the  "terrain  enveloping"  concept  in  which  topog¬ 
raphy  is  calculated  at  a  high  horizontal  resolution  to 
incorporate  the  effects  of  the  sub-grid  scale  features  into 
the  topographic  field.  A  thorough  discussion  of  NORAPS  is 
provided  by  Hodur  (1982  and  1984). 

The  four  major  components  of  NORAPS  are  the  analysis, 
initialization,  forecast  and  output.  The  analysis  component 
consists  of  acquiring  different  types  of  data  (radiosonde, 
pibal,  aircraft,  satellite,  land  and  ship  reports)  and 
applying  quality  control  checks  to  determine  data  validity. 
A  single  bad  observation  can  have  an  adverse  effect  on  the 
regional  model  if  not  removed  prior  to  initialization.  The 
next  step  in  the  analysis  is  to  interpolate  the  observations 
to  the  grid.  The  data  fields  for  the  model  are  the  u  and  v 
wind  components,  temperature,  geopotential,  specific 
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humidity,  surface  pressure,  sea-surface  temperature  and 
terrain  height.  A  regional  update  cycle  is  used  where  the 
12  hour  NORAPS  forecast,  which  serves  as  the  first  guess  and 
the  latest  observations  are  blended  together.  A  sucessive 
corrections  technique  is  then  used  for  the  NORAPS  objective 
analysis  to  improve  the  first-guess  fields  of  wind,  tempera¬ 
ture,  and  geopotential.  An  exponential  weighting  function 
takes  into  account  the  distance  of  the  observation  to  the 
grid  point  and  the  number  of  observations  surrounding  a 
certain  grid  point.  The  analyses  are  performed  at  1000, 
925,  850,  700,  500,  400,  300,  250,  200,  150  and  100  mb. 

The  purpose  of  the  initialization  phase  is  to  approxi¬ 
mately  balance  the  mass  and  wind  fields  and  thus  suppress 
the  growth  of  large  amplitude  inertial-gravity  waves  which 
would  contaminate  the  forecast  fields.  The  static  initiali¬ 
zation  procedure  uses  diagnostic  constraints  to  relate  the 
wind  and  mass  fields.  The  nondivergent  (rotational)  wind 
component  is  obtained  from  the  mass  field  through  the 
balance  equation,  while  the  divergent  (  irrotational )  compo¬ 
nent  is  computed  from  the  omega  equation. 

The  forecast  component  is  the  heart  of  the  NORAPS  model 
and  this  phase  requires  the  majority  of  the  computation 
time.  The  model  uses  the  flux  form  of  the  primitive  equa¬ 
tions  on  a  staggered  grid  scheme  C  (Arakawa,  1977).  This 
grid  scheme  has  excellent  geostrophic  adjustment  properties 
and  group  velocity  characteristics  (Haltiner  and  Williams, 
1980).  The  vertical  coordinate  is  sigma,  which  orients  all 
coordinate  surfaces  parallel  to  the  terrain  surface.  Thus 
the  vertical  velocity  is  identically  zero  at  the  lower 
boundary,  even  in  the  vicinity  of  mountainous  terrain.  The 
vertical  structure  of  the  atmosphere  is  normally  represented 
in  12  discrete  layers,  although  as  many  as  19  layers  may  be 
specified.  All  prognostic  variables  (u,  v,  q  and  T)  except 
vertical  velocity  are  carried  at  the  middle  of  each  layer. 
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A  split-explicit  time  integration  scheme  is  used  to  permit 
larger  time  steps  for  the  slower  meteorological  modes,  while 
still  being  able  to  predict  all  the  gravity  modes.  The  size 
of  the  time  step  is  governed  by  the  computational  stability 
criterion  for  the  horizontal  resolution  selected.  The 
momentum,  thermodynamic  and  moisture  equations  are  solved 
with  conventional  centered  time  ( leapfrog)  and  space  differ¬ 
encing  scheme.  Fourth-order  advection  is  used  for  the 
predictions  equation  set  to  reduce  errors  in  phase  speed. 

One-way  influence  boundary  conditions  are  used  to 
specify  the  time-dependent  lateral  boundary  conditions  on 
the  finer  mesh  NORAPS  model  from  the  NOGAPS  predictions. 
The  one-way  influence  refers  to  the  NOGAPS  solution  forcing 
the  fine-mesh  model,  without  the  fine  grid  affecting  the 
coarse  grid  solution.  For  timeliness  required  in  opera¬ 
tional  use,  these  boundary  conditions  must  be  derived  from 
an  earlier  forecast  rather  than  utilizing  the  corresponding 
NOGAPS  forecast  from  the  same  time.  A  method  developed  by 
Perkey  and  Kreitzberg  (1976)  is  used  to  spatially  interpo¬ 
late  the  solutions  near  the  boundary  of  the  finer  mesh.  The 
NOGAPS  time  tendencies  are  blended  with  the  NORAPS  time 
tendencies  over  a  distance  of  several  grid  points  to  dampen 
spurious  reflections  at  the  regional  model  boundary  due  to 
the  change  in  grid  spacing. 

The  output  phase  of  NORAPS  prepares  the  forecast  data 
for  interpolation  to  the  standard  pressure  levels.  The 
output  fields  can  include  winds,  temperature,  specific 
humidity,  surface  pressure,  relative  humidity,  absolute 
vorticity,  divergence,  surface  sensible  and  latent  heat 
flux,  terrain  height  and  precipitation.  The  output  domain 
can  be  either  the  entire  grid  or  a  subset  of  it.  Three 
different  projections  (Mercator,  Lambert  conformal,  or  polar 
stereographic)  are  available  to  minimize  distortion  in  the 
tropics,  mid-latitudes  or  polar  regions. 


2.  MODEL  PHYSICS 

The  model  physics  contained  in  NORAPS  constitute  a 
crucial  component  in  this  experiment.  The  treatment  of 
diabatic  processes  is  important  to  simulate  the  effect  of 
the  surface  fluxes  across  the  air-sea  interface  on  the 
atmosphere.  NORAPS  includes  representations  of  the 
following  physical  processes: 

(i)  dry  convective  adjustment; 

(ii)  surface  heat,  momentum  and  moisture  fluxes; 

(iii)  cumulus  parameterization; 

(iv)  large-scale  precipitation;  and 

(v)  radiative  transfer  processes 


The  planetary  boundary  layer  (PBL),  which  is  defined 
as  the  lowest  layer  in  the  model  atmosphere,  is  well  mixed 
in  temperature,  momentum  and  moisture.  Interactions  occur¬ 
ring  between  the  lower  boundary  and  overlying  atmospheric 
layer  provide  sources  and  sinks  for  momentum,  heat  and  mois¬ 
ture.  The  effects  of  the  PBL  should  be  included  in  any 
numerical  model  to  physically  simulate  maritime  cyclogenesis 
on  the  time  scales  of  more  than  a  few  hours  as  stated  by 
Anthes  et.  al.  (  1983). 

The  NORAPS  PBL  parameterization  follows  Deardorff 
(  1972).  After  the  layer  mean  values  of  V,  0  and  q  are 
known,  a  bulk  Richardson  number  ( Ri^  in  A. 1)  is  computed  to 
determine  the  stability  of  the  PBL, 


Rih  =  gh(0v  -  0V)/U  2  , 


(A.  1) 


where  g  is  the  gravitational  constant,  h  is  the  PBL  height, 
Qv  is  the  virtual  potential  temperature,  subscript  s  denotes 
surface  values  and  subscript  m  denotes  mean  PBL  values.  For 
unstable  conditions  (Ri]_)<0),  i.  e.  ,  strong  winds,  daytime 
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heating  over  land  surfaces  and  strong  mixing,  a  predictive 
equation  for  the  PBL  height  (h)  proposed  by  Stull  (1976)  is 
used.  The  rate  of  change  of  the  PBL  height  is  related  to 
the  surface  sensible  heat  flux,  mean  PBL  wind  speed,  the 
large-scale  vertical  motion  and  cloud-induced  subsidence. 
For  stable  or  neutral  conditions  (Rib>0),  i.e. ,  light  winds, 
nighttime  over  land  with  weak  mixing,  a  predictive  equation 
for  the  PBL  height  after  Nieuwstadt  and  Tennekes  (1981)  is 
used.  Transfer  (drag)  coefficients  for  heat  (Cg)  and  fric¬ 
tion  (Cu)  are  computed  from  empirical  formulas  that  include 
the  stability  dependence  via  the  Ri^.  Surface  fluxes  of 
heat,  moisture  and  momentum  (A. 2,  A. 3 ,  and  A. 4)  are  computed 
using  the  bulk  aerodynamic  formulas,  which  assume  the 

transfer  coefficients  are  constant  in  time: 

Sensible  heat  flux, 

HS  =  pCp(w'fl')  =  pcpu*cg(0s-0m)  ;  (A.  2) 

Latent  heat  flux, 

H1  =  lv(w'<3')  =  pcpu*cg(qs-qm);  and  (  A.  3 ) 

Surface  stress, 

t  =  pu*2  =  ?um2cu2  ,  (A. 4) 

where  p  is  the  density  of  air,  u*  is  the  frictional 

velocity,  0  is  the  potential  temperature,  Lv  is  the  latent 
heat  of  vaporization,  w  is  the  vertical  velocity,  q  is  the 
specific  humidity,  cp  is  the  specific  heat  of  dry  air  at 

constant  pressure,  subscript  m  denotes  the  mean  PBL  value 
and  subscript  s  denotes  the  surface  value.  The  PBL  is 

constrained  to  extend  through  at  least  the  bottom  layer  of 
the  model  ( approximately  40  mb)  to  avoid  extrapolation  prob¬ 
lems  in  determining  mean  PBL  quantities. 


Another  feature  of  the  NORAPS  PBL  is  that  seasonally 
dependent  climatological  values  of  albedo,  sea  ice,  ground 
wetness  and  surface  roughness  are  specified.  A  predictive 
equation  for  the  ground  temperature  after  Blackadar  (1979) 
is  used  to  model  the  lower  boundary  condition  for  the 
temperature  over  land.  The  sea-surface  temperatures  are 
assumed  to  be  constant  over  the  forecast  period,  which  is 
valid  for  short-range  forecasts. 

2.  Cumulus  Parameterization 

The  NORAPS  model  uses  a  modified  version  of  the  Kuo 
(1965)  cumulus  parameterization  scheme.  This  version  links 
the  convection  to  the  PBL  by  requiring  moisture  convergence 
in  the  PBL.  By  constrast,  the  original  Kuo  version  required 
net  moisture  convergence  in  the  entire  column  before  convec¬ 
tion  was  initiated.  The  moisture  convergence  ( A. 5 )  is 


Mt  =  l/g  V  •  <W<l-ffpbl)  ♦  ps(w'q')s  -  <A-5> 

where  the  first  term  on  the  RHS  is  the  vertically-integrated 
moisture  convergence  and  the  second  term  is  the  surface 
moisture  flux.  Convection  is  assumed  to  occur  when  Mt  >  4.  0 
x  10-D  gm  m-^  s  and  the  equivalent  potential  temperature 
decreases  with  height  (conditionally  unstable)  from  the  PBL 
to  the  first  model  layer  above  the  PBL.  The  final 
constraint  is  that  deep  convection  cannot  occur  if  the 
lifting  condensation  level  (LCL)  is  above  the  PBL.  This 
scheme  partitions  the  moisture  transport  into  two  fractions: 
the  first  (bMt)  serves  to  moisten  the  environment  to  satura¬ 
tion  conditions  through  the  cloud  layer,  and  the  second 
((l-b)Mt)  condenses  and  falls  instantaneously  as  rain.  The 
factor  b  is  expressed  as  the  vertical  average  of  one  minus 
the  relative  humidity  (RH),  where  0  <  RH  <  1.  This  factor 
serves  to  moisten  the  column  when  RH  is  low  and  condense 
moisture  when  the  RH  is  high. 


In  addition  to  large-scale  advection,  temperature 
and  moisture  changes  at  any  level  are  caused  by  convective 
clouds.  It  is  assumed  that  the  temperature  of  the  cloud  is 
warmer  than  the  environment.  Cloud  production,  which  is  the 
ratio  of  the  water  vapor  available  to  the  water  vapor  needed 
to  form  the  cloud,  is  computed  for  each  gridpoint.  The 
fractional  cloud  area  is  used  to  adjust  the  layer  mean 
temperatures  and  moistures  to  account  for  the  presence  of 
clouds. 

3.  Precipitation 

Large-scale  precipitation  ( non-convective )  can  occur 
when  supersaturation  is  achieved  at  any  level.  The  excess 
moisture  is  allowed  to  fall  into  the  next  layer  and  increase 
the  moisture  content  of  that  layer,  or  continue  to  fall  if 
that  layer  is  already  supersaturated.  Precipitation  occurs 
only  when  the  air  is  saturated  from  the  cloud  to  the  ground. 
Convective  precipitation  occurs  according  to  the  modified 
Kuo  cumulus  convection  scheme  discussed  in  the  previous 
section.  The  precipitation  routines  are  only  called  every 
eight  time  steps  for  computational  efficiency.  The  heating 
and  moistening  rates  are  then  spread  evenly  over  subsequent 
time  steps  until  the  next  call  to  these  routines. 

4.  Radiatipn 

The  incorporation  of  solar  radiation  into  numerical 
models  is  essential  for  prediction  of  surface  temperatures 
and  the  cooling  rates  at  cloud  tops  that  may  deepen  cloud 
layers.  The  radiation  parameterization  in  NORAPS  follows 
Katayama  (1974)  for  short-wave  radiation  and  Sasamori  (1968) 
for  long-wave  radiation. 


DATA  ACQUISITION  AND  PROCESSING 

1.  DATA  ACQUISITION 

The  NORAPS  analyses  and  forecast  fields  for  this  budget 
study  were  obtained  from  Dr.  Rich  Hodur  at  the  Naval 
Environmental  Prediction  and  Research  Facility  (NEPRF).  The 
data  on  9-track  tape  were  transferred  to  the  mass  storage 
device  on  the  NPS  IBM  3033.  The  unprocessed  NORAPS  data 
fields  were  on  a  109  x  82  grid  at  12  sigma  levels.  A  hori¬ 
zontal  grid  spacing  of  80  km  is  used.  A  slightly  smaller 
window  (  103  x  76  with  11  pressure  levels)  was  extracted  to 
accomodate  easier  storage  and  access  from  the  disk.  The 
analysis  base  time  for  the  model  run  is  00  GMT  28  March 
1984.  A  36-hour  model  forecast  is  produced  with  output 
fields  generated  every  six  hours  to  12  GMT  29  March  1984. 
NORAPS  analyses  were  available  every  12  hours  from  00  GMT  28 
March  until  12  GMT  29  March  1984. 

The  Lambert  conformal  projection  used  for  the  output 
fields  is  ideally  suited  for  mid-latitudes  since  there  is 
minimum  distortion  between  the  true  parallels  of  30°  and 
60°N.  The  Lambert  conformal  map  is  a  bi-conic,  secant  type 
of  projection  which  preserves  angles  when  projecting  the 
earth's  surface  onto  a  plane  surface. 

To  obtain  the  data  for  the  budget  programs,  several 
preliminary  steps  were  necessary.  First,  surface  and  upper 
air  fields  were  plotted  using  DISSPLA,  which  is  a  software 
package  available  on  the  NPS  IBM  3033  mainframe.  The  user 
must  be  aware  that  DISSPLA  requires  a  rectangular  region  of 
latitude/longitude  points  to  be  specified.  Because  NORAPS 
output  fields  are  specified  in  Lambert  conformal  coordi¬ 
nates,  DISSPLA  would  perform  what  amounts  to  a  double 


transformation  and  produce  a  distorted  and  inaccurate  field. 
Three  separate  steps  are  required  in  the  program  DISPLA 
NORAPS  to  produce  a  plot  with  correct  positioning  on  the  map 
projection.  First,  a  subplot  area  is  specified  and  a 
blanking  routine  is  used  to  truncate  the  lower  curved 
boundary  which  is  standard  for  the  Lambert  conformal  plot. 
This  step  merely  serves  to  ensure  a  rectangular  plot  is 
produced.  Second,  the  contouring  and  a  border  are  drawn 
that  is  separate  from  the  projection  and  geography  routines. 
Third,  the  Lambert  conformal  projection,  i.  e.  ,  the  latitiude 
and  longitude  lines,  and  the  geography  are  added.  Integral 
to  each  of  these  three  separate  steps  is  the  statement,  CALL 
ENDGR(O),  which  terminates  that  particular  block  of  code. 
This  statement  ends  a  subplot  but  remains  on  the  same  phys¬ 
ical  page,  which  allows  other  plots,  such  as  the  contouring 
and  projection  in  this  case,  to  be  drawn  on  the  same  phys¬ 
ical  page. 

It  is  sometimes  necessary  to  retrieve  data  from  the 
output  grid  to  determine  the  center  of  a  low  center  in  the 
(i,j)  Lambert  conformal  coordinates.  The  NDATA  FORTRAN 
program  can  be  easily  modified  to  retrieve  output  data  for 
any  user  specified  field  and  level.  Once  the  appropriate 
data  fields  are  obtained,  a  program  called  TRANS  FORTRAN  is 
used  to  transform  the  (i,j)  Lambert  conformal  coordinate  to 
a  latitude  and  longitude  on  the  earth's  surface.  The  corre¬ 
sponding  latitude  and  longitude  of  the  low  center  are 
entered  into  an  interactive  program  (called  STORMO  FORTRAN) 
at  each  time  period  to  compute  the  speed  and  direction  of 
the  cyclone  center.  A  forward  difference  is  used  to  compute 
the  speed  for  the  first  and  last  time  periods,  while  a 
centered  difference  is  used  for  the  other  time  periods.  The 
low  center  location,  direction  and  speed  are  then  entered  at 
the  end  of  the  budget  programs  as  required  parameters  for 
the  budget  calculations. 


2.  DESCRIPTION  OF  WIND  ADJUSTMENT 


Because  all  NORAPS  output  fields  are  given  in  Lambert 
conformal  coordinates,  the  wind  directions  will  be  distorted 
from  the  true  direction  away  from  the  central  meridian  of 
the  conformal  grid.  This  central  or  true  meridian,  which  is 
80°W  in  this  case,  is  parallel  to  the  y-axis  in  an  x,y 
Cartesian  coordinate  system.  A  subroutine  called  WNDADJ 
makes  the  necessary  transformations  to  provide  the  true  wind 
direction  on  the  earth's  surface.  The  only  information 
required  to  make  this  transformation  are  the  Lambert 
conformal  coordinates  of  the  pole  point,  which  are  52. 0 
(  x-coordinate )  and  -37.986  {  y-coordinate ) .  Equations  B.  1 

and  B.  2  are  used  to  transform  the  Lambert  conformal  wind 
direction  components,  denoted  u'  and  v'  are: 

U  =  u'cos(0)  -  v'sin(0);  and  (  B.  1) 


V  =  v'cos(0)  +  u'sin(0)  ,  (  B.  2 ) 

where  U  and  V  are  the  true  horizontal  wind  direction  compo¬ 
nents  and  the  angle  0  =  tan_1(xp-  x/y  -  yp),  where  xp,  yp 
are  the  Lambert  conformal  pole  point  coordinates  and  x, y  are 
the  Lambert  conformal  (i,j)  grid  point  that  is  being  trans¬ 
formed.  The  convention  for  the  angle  (0),  which  the  Lambert 
conformal  coordinate  axes  are  rotated  to  become  true  direc¬ 
tional  axes,  is  positive  for  a  counter-clockwise  direction 
and  negative  for  a  clockwise  rotation. 

3.  CONVENTIONS  FOR  NORAPS  FIELDS  AND  BUDGET  PROGRAMS 

The  NORAPS  output  fields  are  presented  in  a  right-hand 
coordinate  system  with  the  (1,1)  grid  point  being  the  south¬ 
west  corner  of  the  grid.  The  column  value  increases  east¬ 
ward  and  the  row  value  increases  northward.  These  NORAPS 


data  fields,  however,  are  read  in  for  the  budget  programs  in 
a  different  manner.  A  left-hand  coordinate  system  is  used 
in  the  budget  program  and  therefore,  data  fields  have  to  be 
read  in  accordingly.  The  fields  are  read  in  using  the  HOJO 
subroutine  with  the  (1,1)  grid  point  at  the  northwest 
corner.  The  column  value  increases  eastward  and  the  row 
value  increases  downward.  The  pole  point  is  referenced  in 
this  coordinate  system  since  wind  adjustments  are  performed 
after  the  data  fields  have  been  read  in. 

The  latitude  convention  in  the  budget  program  is  posi¬ 
tive  north  and  negative  south;  longitude  is  positive  west  of 
Greenwich  and  (360-longitude),  east  of  Greenwich.  The 
convention  for  the  normal  wind  components  in  the  budget 
programs  are  positive  outward  and  negative  inward.  This 
likewise  applies  to  the  normal  component  of  the  cyclone 
velocity. 

Another  convention  the  user  should  be  aware  of  is  the 
method  of  defining  latitude  and  longitude  in  the  DISSPLA 
software  package  used  in  plotting  the  NORAPS  output  fields. 
Longitudes  west  of  Greenwich  are  negative,  while  east  of 
Greenwich  they  are  positive.  The  latitudes  are  positive  for 
north  and  negative  for  south.  These  are  arguments  to  be 
included  in  the  subroutines  GRAF  and  MAPGR. 


NORAPS  MODEL  RESULTS  -  CORRECTED  VERSION 


1.  GENERAL 

The  incorrect  specification  of  moistening  and  heating 
rates  in  the  NORAPS  model  led  to  a  poor  comparison  of  budget 
results  with  the  model-predicted  diabatic  heating  and  mois¬ 
tening  rates  as  discussed  in  Chapter  5.  The  corrections  to 
the  specification  of  the  heating  and  moistening  rates  in  the 
prognostic  equations  were  made  in  a  test  version  of  the 
NORAPS  model  and  the  forecast  was  re-run  for  the  28-29  March 
1984  storm.  Corrections  to  the  calculation  of  the  vertical 
fluxes  of  moisture  and  momentum  were  also  made  so  that  a 
linear  (in  pressure)  interpolation  scheme  is  used  to  esti¬ 
mate  horizontal  winds  and  specific  humidity  at  levels  of 
vertical  velocity.  An  abbreviated  discussion  of  the  results 
in  the  corrected  model  run  is  presented.  A  brief  synopsis 
is  given  for  the  track  and  intensity  of  the  low  center  in 
the  two  versions  of  the  model.  The  new  residuals  for  the 
heat  and  moisture  budgets  are  then  compared  with  the  updated 
model-predicted  heating  and  moistening  rates.  These  model- 
predicted  heating  and  moistening  rates  are  accumulated  6  h 
values  converted  to  a  daily  rate  rather  than  an  instanta¬ 
neous  value  obtained  each  6  h  during  the  integration  of  the 
model  and  converted  to  a  daily  rate.  Results  are  presented 
for  radius  4°  lat.  as  before  and  figures  are  placed  at  the 
end  of  this  appendix. 

2.  SYNOPTIC  DISCUSSION 

The  track  and  intensity  of  the  surface  low  center  during 
the  first  24  h  in  the  corrected  NORAPS  model  run  are  very 
similar  to  the  incorrect  version  (Fig.  C.  la).  A  weaker 
intensity  low  center  is  predicted  in  the  corrected  model 


during  the  last  12  h.  A  slightly  greater  cyclone  speed  is 
predicted  in  the  revised  model  during  the  last  12  h,  which 
reduces  the  position  error  at  12  GMT  29  March  1984  (Fig. 
C.  lb). 

During  00-06  GMT  29  March  1984,  the  sea-level  pressure 
(SLP)  tendency  diminishes  from  4.5  mb/6  h  (old  version)  to 
0.9  mb/6  h  (new  version).  During  06-12  GMT  29  March  1984, 
the  SLP  tendency  only  slightly  diminishes  from  3.  7  mb/6  h  to 


3.0  mb/6  h.  Thus,  the  central  SLP  in  the  new  version  is  978 
mb  at  12  GMT  29  March  1984,  whereas  the  intensity  of  the  low 
center  is  974  mb  in  the  old  version.  By  comparison,  the 
NORAPS  analysis  at  12  GMT  29  March  1984  has  a  973  mb  center 
after  a  8.  5  mb/12  h  SLP  decrease,  which  is  nearly  double  the 
SLP  tendency  for  the  new  version.  Although  a  weaker  storm 
results  from  the  incorporation  of  corrections  to  the  NORAPS 
model,  the  position  error  improves  at  36  h. 

A  noteworthy  difference  between  the  old  and  new  forecast 
is  the  formation  of  a  secondary  low  pressure  center  within 
the  northeast  quadrant  of  the  major  center  of  low  pressure 
at  00  GMT  29  March  1984.  This  mesoscale  feature  has  some 
similar  characteristics  as  the  coastal  low  that  becomes  the 
major  low  center  during  the  last  12  h  of  the  forecast 
period.  The  exact  details  of  this  mesoscale  feature  are  not 
correctly  predicted  by  the  N0RAPS  model,  but  it  presents  a 
much  more  realistic  picture  of  the  actual  situation  as  seen 
in  the  NMC  analyses  discussed  in  Chapter  3. 


3. 


BUDGET  RESULTS 
1.  Heat  Budget 

The  terms  in  the  heat  budget  for  the  new  forecast 
have  a  similar  vertical  profile  as  in  the  old  forecast.  A 
noticeable  difference  in  the  new  forecast  is  the  greater 
oscillatory  trend  in  time  of  the  vertical  temperature  advec- 
tion  and  energy  conversion  terms  (not  shown).  This  is  a 
result  of  a  stronger  secondary  maximum  in  the  vertical 
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velocity  field  (Fig.  C.  2 )  at  09  GMT  28  March  1984  in  the  new 
forecast. 

The  residual  (diabatic  heating  plus  calculation 
error)  diagnosed  from  the  budget  for  the  new  forecast  ( Fig. 
C. 3a)  has  two  maxima:  a  20°C/day  center  at  450  mb  between 
09  and  15  GMT  28  March  1984  and  a  25°C/day  center  at  09  GMT 
29  March  1984  near  350  mb.  The  magnitude  of  maximum  heating 
is  5°C/day  greater  than  in  the  old  forecast  and  the  level  is 
100-200  mb  higher.  The  convective  precipitation  component 
does  not  increase  substantially  during  the  later  stages  of 
the  storm  and  therefore  can  not  account  for  this  higher 
level  of  heating.  It  had  been  hoped  that  the  diabatic 
heating  would  be  at  lower  levels  using  the  new  forecast 
fields.  The  diagnosed  heating  rate  (Fig.  C.3a)  in  the  new 
forecast  is  much  closer  to  the  actual  model  heating  rate 
(Fig.  C. 3b).  The  magnitude  of  20-25°C/day  during  the  first 
strong  heating  period  at  09  GMT  28  March  1984  agrees  well 
with  the  actual  rate  in  the  model.  Excellent  agreement  is 
also  found  in  the  phases  and  levels  of  maxima  in  the  diag¬ 
nosed  and  predicted  heating  rates.  Stronger  heating  occurs 
at  the  last  time  period  than  is  diagnosed  in  the  heat 
budget.  Other  differences  between  the  diagnosed  and 
predicted  rates  are  the  higher  level  of  heating  (350  mb) 
diagnosed  in  the  heat  budget  at  09  GMT  29  March  1984  and  the 
strong  heating  (  in  the  budget)  below  850  mb  between  09  and 
21  GMT  28  March  1984  as  opposed  to  cooling  predicted  by  the 
model  during  this  period. 

The  difference  between  the  budget  diagnosed  heating 
rate  and  the  model-predicted  rate  is  given  in  Fig.  C.  4. 
Errors  of  less  than  20-25%  in  the  budget  are  acceptable  and 
can  be  attributed  to  computational  errors.  If  the  errors 
are  larger  than  this,  then  the  validity  of  the  particular 
budget  must  be  questioned.  Two  regions  have  significant 
errors  (  greater  than  10°C/day)  in  the  diagnosed  heating 
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rate.  The  first  region  occurs  below  850  mb  centered  at  15 
GMT  28  March  1984.  The  large  surface  heat  flux  ( 10°C/day) 
at  this  time  can  explain  this  difference  since  this  surface 
heat  flux  has  not  been  included  in  the  model-predicted 
diabatic  heating  rate.  The  second  region  of  large  errors 
(negative)  occurs  at  the  last  time  period  in  the  middle 
troposphere.  These  errors  are  probably  due  to  the  fact  that 
the  largest  heating  rates  are  predicted  by  the  model  at  this 
time  and  the  time-averaging  process  in  the  budget  tends  to 
weaken  the  heating  signal  by  spreading  the  maximum  across 
two  time  periods.  An  alternate  explanation  is  the  strong 
vertical  velocity  diagnosed  at  the  last  time  period  is  asso¬ 
ciated  with  excessive  adiabatic  cooling,  which  is  not 
balanced  by  the  other  terms  in  the  budget.  Errors  of 
5-10°C/day  above  350  mb  are  due  to  unrealistically  large 
values  of  horizontal  and  vertical  temperature  advection  at 
the  upper  levels  in  the  diagnostic  calculations. 

2.  Moisture  Budget 

The  moisture  budget  for  the  new  forecast  exhibits  a 
similar  structure  as  in  the  old  version  with  the  exception 
of  the  vertical  moisture  flux  term.  This  term  has  much  less 
vertical  structure  after  21  GMT  28  March  1984  (Fig.  C. 5)  as 
a  result  of  the  corrections  in  the  moistening  rate.  In 
particular,  realistic  values  of  specific  humidity  are  now 
found  in  the  middle  to  upper  troposphere. 

The  budget  residual  for  the  new  forecast  (Fig.  C. 6a) 
has  (negative)  maxima  in  moistening  at  09  GMT  28  March  and 
09  GMT  29  March  1984.  These  maxima  at  700  mb  and  600  mb, 
respectively,  compare  rather  poorly  with  the  levels  of 
maximum  diabatic  heating  (Fig.  C. 3a),  which  are  diagnosed  to 
be  200-300  mb  higher.  In  the  new  forecast,  the  actual  model 
moistening  rate  (Fig.  C. 6b)  exhibits  a  peak  magnitude  in 
moistening  at  450-500  mb.  This  compares  rather  poorly  with 
the  budget  diagnosed  results,  which  are  150-250  mb  lower  in 
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the  troposphere.  However,  excellent  agreement  exists  in  the 
phases  and  maximum  magnitudes  of  the  moistening.  The  lower 
elevation  of  the  residual  in  the  moisture  budget  is 
distinctly  different  from  the  heat  budget  where  the  diag¬ 
nosed  maximum  heating  occurs  at  a  higher  level  than  in  the 
model. 

The  differences  between  the  budget  diagnosed  and 
model-predicted  moistening  rate  are  illustrated  in  Fig.  C. 7. 
Large  errors  (positive  and  negative)  occur  throughout  the 
depth  of  the  troposphere.  The  negative  errors  indicate  an 
overestimation  of  the  moisture  sink  in  the  budget  residual, 
while  positive  values  indicate  an  underestimation  of  the 
source.  These  errors  in  the  moisture  budget  are  on  the 
order  of  the  magnitude  of  the  moistening  rates,  especially 
at  09  GMT  28  March  and  09  GMT  29  March  1984.  Evidently, 
there  is  another  process  in  the  forecast  model  which  redis¬ 
tributes  moisture  in  the  vertical  but  is  not  included  in  the 
moistening  rate  we  have  output  from  the  model.  Other  poten¬ 
tial  sources  of  error  are  the  highly  variable  (horizontal 
and  vertical)  moisture  fields,  too  coarse  of  a  temporal 
resolution  and  unresolved  processes  in  this  rapidly  varying 
situation. 

The  agreement  between  the  column-integrated  budget 
and  model-predicted  moistening  rates  is  significantly 
improved  in  the  new  forecast.  In  the  old  forecast,  a  factor 
of  two  difference  was  found  (Fig.  C. 8a)  between  the  two 
rates.  The  diagnosed  residual  (Fig.  C. 8b)  has  a  non- 
systematic  error  in  the  new  forecast  which  is  in  much  closer 
agreement  with  the  model-predicted  moistening  rate.  Thus, 
the  correction  to  the  prediction  of  the  moistening  rate  in 
the  model  leads  to  consistent  results  between  the  budget 
diagnosed  and  model-predicted  moistening  rates  in  a  column- 
averaged  sense.  Another  pleasing  result  is  the  degree  of 
coincidence  of  the  actual  model-moistening  rate  and  the 


model  precipitation  when  accumulated  values  rather  than 
instantaneous  rates  are  calculated  for  both  variables. 

4.  SUMMARY 

In  general,  the  heat  budget  captures  rather  well  the 
major  thermodynamic  process  of  diabatic  heating  in  compar¬ 
ison  to  the  rates  internally  predicted  within  the  model. 
The  corrections  to  the  heating  and  moistening  rates  in  the 
the  NORAPS  model  have  yielded  more  consistent  results  in  the 
column-integrated  moisture  budget.  However,  a  discrepancy 
remains  in  the  vertical  distribution  of  moistening  estimated 
from  the  budget  and  directly  from  the  model  prediction. 

Problems  still  remain  in  the  model  and  budgets  which 
need  to  be  resolved.  The  most  serious  problem  is  the 
difference  in  the  elevation  of  the  maximum  heating/ 
moistening  rates  between  the  budget  and  model-predicted 
calculation.  Unrealistic  vertical  temperature  profiles 
above  350  mb  are  found  in  the  model  output  data. 
Unrealistically  strong  horizontal  and  vertical  temperature 
advection  above  350  mb  is  diagnosed  in  the  heat  budget.  The 
moisture  budget  may  be  less  reliable  than  the  heat  budget 
due  to  the  higher  degree  of  noise  in  the  specific  humidity 
field  as  compared  to  the  temperature  field. 
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